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ABSTRACT 

The Groth field is one of the sky regions that will be targeted by the OTELO (OSIRIS Tunable 
Filter Emission Line Object) survey in the optical 820 nm and 920 nm atmospheric windows. This 
field has been observed by AEGIS (All-wavelength Extended Groth strip International Survey) 
covering the full spectral range, from X-rays to radio waves. Chandra X-ray data with total 
exposure time of 200ksec are analyzed and combined with optical broadband data of the Groth 
field in order to study a set of structural parameters of the X-ray emitters and its relation with X- 
ray properties. We processed the raw, public X-ray data using the Chandra Interactive Analysis 
of Observations and determined and analyzed different structural parameters in order to produce 
a morphological classification of X-ray sources. Finally, we analyzed the angular clustering of 
these sources using 2-point correlation functions. We present a catalog of 340 X-ray emitters 
with optical counterpart. We obtained the number counts and compared them with AEGIS data. 
Objects have been classified by nuclear type using a diagnostic diagram relating X-ray-to-optical 
ratio (X/O) to hardness ratio (HR). Also, we combined structural parameters with other X-ray 
and optical properties, and found for the first time an anticorrelation between the X/O ratio 
and the Abraham concentration index which might suggest that early type galaxies have lower 
Eddington rates than those of late type galaxies. A significant positive angular clustering was 
obtained from a preliminary analysis of 4 subsamples of the X-ray sources catalog. The clustering 
signal of the opticaly extended counterparts is similar to that of strongly clustered populations 
like red and very red galaxies, suggesting that the environment plays an important role in AGN 
phenomena. 

Subject headings: Galaxies: active — X-rays: galaxies 



Introduction 



The OTELO project ijCena et al.l I200R l2008h 
is a flux-limited survey of emission-line objects 
in large and perfectly defined volumes of the 
Universe. The OTELO survey, more than one 
magnitude deeper than other emission line sur- 
veys, includes a wide variety of astronomical 
sources: star-forming galaxies, starburst galax- 
ies, emission-line ellipticals, AGNs, QSOs, Lya 
emitters, peculiar stars, etc. With these data, a 
number of scientific problems will be addressed, 



including: star formation rates (SFR) and metal- 
licity evolution of galaxies; AGN and QSO evo- 
lution. Deep imaging progra ms such as sev- 
eral VLT surveys for GOODS foickinson et~al 



20041). COSMOS (ICapaket alJ l2007h and SXDS 
( Furusawa et al. 2008f ) fields will serve as prepara- 



tory work for OTELO. Likewise, we have per- 
formed a deep BVRI broadband survey in one 
of the fields of the Groth strip selected for the 
OTELO survey (|Cepa et al.ll2008h . 

In the present work, data from this broad band 
survey is combined with other wavelength data 
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for achieving several objectives of the OTELO 
project. Specifically, the matching of these op- 
tical data with X-ray allows succesfully tackling 
the study of t he AGN population as shown i n 
Barcons et al.l (120071). iGeorgakakis et alj (|2006l ). 
and ISteffen et alj ~( 20061 ). Optical imaging sur- 
veys allow a morphological classification of all 
sufficiently resolved sources. This is esse ntial to 
reveal the nature of AG N host galaxies (e.g. Hasanl 
2007tlFiore et al.ll200l[ ). although the morphologi- 
cal classification becomes difficult at high redshift, 
due to the reduced resolution of images. Several 
work s J Graham et al.l 20011: ; Ferrarese fc Merrill 
200d : IGebhardt etaU l2000h have shown that 
the AGNs are directly related with some host 
galaxy properties , parti cularly with their bulges. 
Kauffmann et al. (2003) have pointed out that 
the AGNs in the local universe are hosted pre- 
dominant ly in bulge-d ominated galaxies. More 
recently, iHasanl (|2007l ) using HST and Chandra 
data of the GOODS South field has found that 
the most moderate luminosity AGNs hosts are 
bulge-dominated in the redshift range z~ 0.4-1.3. 
A tight correlation between black hole mass and 
the bul ge dispersion velocity has been well estab- 
lished (jFerrarese fc Merritl I2OO0I IGebhardt et~al 



2000) , as well as with the concentration of bulges 
([Graham et"aLll2001bh . 

On the other hand, deep X-ray extragalactic 
surveys are extremely useful tools to investigate 
the time evolution of the population of active 
galactic nuclei (AGN) and to shed light on their 
triggering mechanisms. The exceptional effective- 
ness at finding AGN arises largely because X- 
ray selection has reduced absorption bias, minimal 
dilution by host-galaxy starlight and allows effi- 
cient optical spectroscopic follow-ups of high prob- 
ability AGN candidates with faint optical coun- 
terparts. The Chandra and XMM-Newton ob- 
servatories have revolutionized this field of re- 
search, generating the most sensitive X-ray sur- 
veys ever performed. In fact, deep surveys pro- 
duced by these observatories have increased the 
resolved fraction of the cosmic X-ray background 
(CXRB) to about 90% in the 0.5-2 keV range. 
A large part of the detected X-ray sources are 
AGN (< 70%), although X-ray emitters include 
a mix of different types of objects such as galaxy 
clusters and stars. On the other hand, these 
powerful observatories make it possible to access 



the hard (2-10 keV) range, allowing to directly 
probe the AGN activity not contaminated by star 
formation processes. Moreover, hard X-ray sur- 
veys are capable to detect all but the most ab- 
sorbed (Compton thick) sources. Therefore, they 
provide th e most complete a n d unbiased samples 
of A GNs (|Mushotzkvj 12004 iBrandt k Hasingerl 
2005). Number counts relations have now been 



determined (e.g. ICappelluti et al.ll2007l ) although 
with some evidence for field-to-field variations. 
Such variations are expected at some level owing 
to "cosmic statistics" associated with large-scale 
structures and clustering feat ures that have been 
detected in the X-ray sky (iBarger et al. .2003; 



Yang et a l. 2003; Gilli et al.ll2003l 12005). 



In this paper we present the analysis of pub- 
lic, deep (200 ksec) Chandra/ ACIS observations 
of three fields comprising the original Groth- 
Westphal strip (GWS), gathered from the Chan- 
dra Data Archive, and combined with optical 
BVRI data from our broadband survey carried 
out with the 4.2m Willia m-Herschel Teles cope 
(WHT) at La Palma (see ICepa et al.l l2008l for 
a detailed description). We have tried to imple- 
ment an innovative approach to the optical study 
of X-ray emitters by investigating correlations 
between broadband optical and X-ray parame- 
ters, specifically the X-ray-to-optical ratio and 
hardness ratios, and optical structural parameters 
that provide information about the host galaxy 
morphology and populations, like asymmetry and 
concentration indexes and optical colors. 

This paper is structured as follows: in Section 
[21 we describe the observational data, including 
X-ray data processing and source detection, as 
well as optical broadband data and the selection 
of the sample of X-ray sources with optical coun- 
terparts. Section [3] reviews the X-ray properties 
of the galaxy sample, the procedures for comput- 
ing optical structural parameters, as well as the 
morphological analysis carried out. The reliabil- 
ity of different structural parameters used in mor- 
phological classification is also discussed. Section 
|4] presents the nuclear type classification, the re- 
lationship between X-ray and optical structural 
properties, and the comparison with results from 
the literature. The catalog of X-ray emitters with 
optical counterparts and their morphological clas- 
sification is presented in the electronic edition. Fi- 
nally, in Section [5] we present the results of the 
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clustering analysis (based on the angular corre- 
lation function formalism) applied on 4 selected 
subsamples of the X-ray catalog. 

2. Observational data 
2.1. X-ray data 

2.1.1. Observations and data processing 

Chandra has observed three consecutive fields 
centered at the original HST Groth-Westphal 
strip (GWS) using the ACIS-I instrument. All 
datasets have been gathered from the Chandra 
Data Archive (http://asc.harvard.edu/cda/) 
using the Chaser tool. PI of all the retrieved 
Chandra observations is K. Nandra. Total expo- 
sure time in each field is about 200 ksec. The total 
field size of ACIS-I chips 0, 1, 2 & 3 is 16.9' x 16.9' 
(ACIS-S chips 2 & 3 were also used, but their 
FOVs do not overlap with our optical data and 
therefore were not used). The FAINT and VFAINT 
telemetry modes have been used, with telemetry 
saturation rates of 170 and 67 events/s, respec- 
tively. Due to the lack of bright sources within 
the fields, we assume that there are no pile-up ef- 
fects in any observation. Table [1] summarizes the 
main characteristics of the Chandra observations. 

We processed the data Using the Chandra 
Interactive Analysis of Observations (CIAO) 
(http : / / exc . harvard. edu/ ciao/), v4.0 and Cal- 
ibration Data Base (CALDB) v3.4.0. Standard 
reduction procedures have been applied, as de- 
scribed in the CIAO Science Threads to produce 
new level 2 event files. The output level 2 event 
files have been restricted to the 0.5-8 keV range 
to avoid high background spectral regions. Fur- 
thermore, we analyzed the source-free light curves 
in order to define additional good time intervals 
(GTI). To this end, we applied the CIAO detection 
program celldetect with a threshold S/N = 3. 
The detected sources have been removed upon 
creation of the time-binned light curves, that have 
been then used for the computation of the GTIs. 
A ±3cr rejection criterion has been applied to re- 
move high background intervals. After filtering 
for the additional GTIs, the level 2 event files be- 
longing to the same target were co-added. To this 
end, we first improved the absolute astrometry 
(the overall 90% uncertainty circle of Chandra X- 
ray absolute position has a radius of 0.6 arcsec) 



Table 1: Chandra observations log for the three 
pointings of the GWS field. Col. (1): observa- 
tion number; Col. (2,3): RA and DEC (J2000); 
Col. (4): Standard Data Processing (SDP) ver- 
sion; Col. (5): Observation format in the Time 
exposure mode; Col. (6): Net exposure time (Sum 
of Good Time Intervals, GTI) (ksec) 



Obs 


RA 


DEC 


SDP 


Mode 


GTI 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


5853 


14:15:22.5 


+52:08:26.4 


7.6.7.1 


VFAINT 


42.4 


5854 






7.6.7.1 


VFAINT 


50.2 


6222 






7.6.7.1 


VFAINT 


35.2 


6223 






7.6.7.1 


VFAINT 


49 


6366 






7.6.7.1 


VFAINT 


15.96 


7187 






7.6.7.1 


VFAINT 


8 


5851 


14:16:24.5 


+52:20:02.59 


7.6.7.1 


VFAINT 


36 


5852 






7.6.4 


VFAINT 


11.95 


6220 






7.6.7.1 


VFAINT 


36 


6221 






7.6.7.1 


VFAINT 


3.8 


6391 






7.6.7.1 


VFAINT 


9.4 


7169 






7.6.4.1 


VFAINT 


18.2 


7181 






7.6.7.1 


VFAINT 


8 


7188 






7.6.4.1 


VFAINT 


4.2 


7236 






7.6.4 


VFAINT 


18.8 


7237 






7.6.4.1 


VFAINT 


17.2 


7238 






7.6.4 


VFAINT 


9.6 


7239 






7.6.4.1 


VFAINT 


16.2 


3305 


14:17:43.0 


+52:28:25.2 


6.8.0 


FAINT 


27.6 


4357 






6.9.0 


FAINT 


79.8 


4365 






6.9.0 


FAINT 


58.2 



applying corrections for translation, scale and ro- 
tation to the world coordinate system (WCS) of 
each data file, by comparing two sets of source 
lists from the same sky region. Average posi- 
tional accuracy improvement is 6.3% with respect 
to the original astrometry provided by the space- 
craft attitude files. Once we have improved the 
astrometry, the event files have been merged. 

The output merged files have been filtered 
to create event files for several energy bands: 
full (0.5-7 keV), soft (0.5-2 keV), hard (2-7 keV) 
hard2 (2-4.5 keV) and vhard (4-7 keV). Unbinned 
effective exposure maps at a single energy, rep- 
resentative of each band, were created at 2.5 keV 
(full), IkeV (soft), 4keV (hard), 3keV (hard2) 
and 5.5 keV (vhard). Resulting maps have units of 
time x effective area, i.e. s cm 2 counts photon -1 . 



4 



2.1.2. Source detection 

We applied the CIAO wavdetect Mexican-Hat 
wavelet source detectiorfj] program to all bands 
(full, soft, hard, hard2 and vhard) for the three 
fields. Several wavelet scales have been applied: 
1, V2, 2, 2\/2, 4, 4 v / 2, 8, 8^2 and 16 pixels. 
Small scales are well suited to the detection of 
small sources while larger scales are appropriate 
for more extended sources. Due to the "blank 
field" nature of the GWS, we do not expect to 
detect sources at a scale larger than 16 pixels (~ 
8 arcsec). For each applied scale, wtransf orm pro- 
duces a complete set of detections. We set a signif- 
icance threshold of 2 x 10 -7 . Since this wavdetect 
input parameter measures the number of spurious 
events per pixel, we expect «0.2 fake detections 
per detector of 1024 x 1024 pixels i.e. about 12 
spurious sources in the entire catalog (4 detectors 
x 5 bands x 3 fields). The output wavdetect 
detection files are fits tables containing the posi- 
tion of sources (RA & DEC), count rates in pho- 
tons s" 1 cm~ 2 and ancillary information. The final 
X-ray source catalog has been obtained by cross- 
matching detections in all five bands in order to 
increase the reliability of detections. Best match is 
given by a minimum separation of sources around 
a great circle. We set a maximum search distance 
of 2 arcsec (see section I2.4|) . We computed the 
hardness ratios as follows: 



HR(A 1 E,A 2 E) = 



CRjAjE) -CR(A 2 E) 
CR{A 1 E) + CR(A 2 E) 



(1) 



where A\E y A 2 E are different energy bands and 
CR (A n E) is the count rate in a given energy 
band. Four hardness ratios have been set in this 
way: HRi = HR(hard, soft), HR'j = HR(hard2, 
soft), HR 2 = HR(very hard, hard) and HR' 2 = 
HR(very hard, hard2). 

We also computed the energetic fluxes from count 
rates. To this end, mean photon energies in 
each band have been calculated assuming a power 
law spectrum with T = 1.5 (a reasonable assump- 
tion after inspecting the hardness ratios diagram 
and diagnosis grid depicted in Figured]), redshift 
z = 0.5, galactic absorption n# = 1.3 x 10 20 cm~ 2 



and no intrinsic absorption. We applied the ab- 
sorption coefficients provided by the TRANNM func- 
tion within the package PIMMS v3.9b (Mukai, 
1993). The output catalog has 639 unique X-ray 
emitters (after removing 20 common sources in the 
field overlapping area). From them, 76 objects are 
located outside the field covered by optical data. 
From the total sources, 46% (297) have been de- 
tected with a significance > 4cr, 21% (132) above 
the 3cr significance level and 13% (81) down to 
2.5cr level. 

The number of sources detected in each band 
are 490 (full), 465 (soft), 284 (hard), 266 (hard2) 
and 121 (vhard), with median fluxes of 1.81 x 10 -15 , 
5.57 x 10~ 16 , 2.34 xlO" 15 , 1.58 x 10~ 15 and 
3.75 x 10 -15 ergcm~ 2 s _1 , and with limiting fluxes 
at 3ct level of 4.8 x 10~ 16 , 1.1 x 10~ 16 , 2.8 x 10~ 16 , 
1.3xl0~ 16 and 7.3 x 10" 16 ergem" 2 s" 1 , respec- 
tively. Median errors in fluxes (for objects above 
the 3ct significance level) are 8% in full band, 10% 
in soft band, 12% in hard band, 14% in hard2 band 
and 22% in vhard band. There are 16 sources 
detected only in hard band, 18 in hard2 band 
and 3 sources are only observed in vhard band. 
We find an overabundance of soft X-ray sources 
in our final catalog. There are 103 sources de- 
tected only in this band. This effect can be at 
least partially explained considering the decrease 
of telescope + ACIS effective area from some 600 
cm 2 in the 1-2 keV range (including the 0.5-2 keV 
band defined as soft), until approximately 200 cm 2 
at ~ 6 keV. On the other hand, the total back- 
ground between 0.5 and 7keV is around 4 times 
higher than in the 0.5-2 keV range. The combina- 
tion of these two factors justifies the high Chan- 
dra/ ACIS source detection efficiency in the soft 
X-ray band and the mentioned overabundance. 
The first version of this cat alog was presented in 



Sanchez-Portal et all (2007) 



The fact that the PSF of X-ray detectors often has 
Gaussian-like shape motivates our use of the Marr wavelet, 
or "Mexican Hat" function. 



2.2. Comparison with AEGIS data 

The All-wavelength Extended Groth strip In- 
ternational Survey (AEGIS) collaboration has re- 
cently made public their data sets, including a 
catalog of X-ray sources that comprises the ob- 
servations used within this paper. The complete 
AEGIS X-ray catalog (|Laird et all 120081 ) consists 
of 1325 sources; from these, 471 are located within 
the area covered by our data (fields EGS-6, 7 and 
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8 according to their nomenclature) . From our 429 
sources detected above 3a significance level, 416 
coincide with AEGIS sources using a search ra- 
dius of 2 arcsec (388 sources if only unique matches 
are considered). Figure [1] shows a good agreement 
between the source counts used in this work and 
those provided by the AEGIS team: median rela- 
tive deviations of our data with respect to AEGIS 
are 2% (full band), 0.7% (soft), 5.5% (hard) and 
7.8% (vhard). 




200 400 



600 1000 1200 

countSfondhis work) 




(counts^nlthis workJ-counts^nfAEGIS)) x lOO/countefunOhis work) 



Fig. 1. — Upper pannel: full band net counts from 
this work vs AEGIS net counts. The solid line rep- 
resents the best linear fit, with a slope of 0.99 and 
zero point 5.23 counts. Lower pannel: histogram 
of relative deviations of our full band count data 
with respect to AEGIS. Median deviation is lower 
than 2%. 



2.3. Optical data 

We observed three pointings in the direction of 
the GWS, each of them using the B, V, R, and I 
broadband filters. The total area covered is 0.18 
square degrees. The observations have been car- 



ried out using the Prime Focus Imaging Platform 
(PFIP) at the 4.2m William-Herschel Telescope 
(WHT) in the Roque de los Muchachos Observa- 
tory (La Palma, Canary Islands). This camera has 
two CCD detectors (2k x 4k), with a total field of 
view of 16' x 16' and a pixel size of 0.237 arcsec. 
Several exposures have been performed at each 
pointing, of 600, 800, 900, or 1000 sec, depending 
on the filter, with a dithering of 15 arcsec between 
consecutive exposures in order to allow eliminat- 
ing cosmic rays and to fill the gap between de- 
tectors. Standard reduction procedures have been 
applied. Resulting limiting magnitudes are 25, 25, 
24.5 and 23.5 in B, V, R and I band, respectively. 
A detailed description of data reduct ion and op- 
tical s ource detection can be found in lCepa et al.l 
(|2008l ) . The optical catalog used in this work con- 
tains data for ~ 44000 objects. 

2.4. The catalog of optical counterparts 

In order to build a sample of X-ray emitters 
with optical counterparts, we performed cross- 
matching of the X-ray and opt ical catalogs using 
the TOPCAT tool (|Tavlorl 120051 ). applying a crite- 
rion of minimum sky distance along a great circle. 
X-ray sources located inside the field covered by 
optical data and their optical counterparts are rep- 
resented in Figure [2 

We chose a maximum search radius of 2 arcsec af- 
ter carrying out several tests considering a maxi- 
mum searching distance of 1 to 5 arcsec with a bin 
of 0.5 arcsec (Figure [3]). As expected, the num- 
ber of counterparts increases with the maximum 
search radius, but also the number of multiple 
matches within the error circle. Setting a search 
radius of 1 arcsec we minimise the number of mul- 
tiple cross-matches, but at the cost of loosing an 
outstanding number of potential counterparts. On 
the other hand, choosing a search radius greater 
than 2 arcsec produces a large number of multiple 
matches, increasing the uncertainty in counterpart 
identification. Therefore, we set a 2 arsec search 
radius as an adequate compromise. This practi- 
cal criterion is consistent with the more rigorous 
analysis performed below. 

We detected 340 optical counterparts of the X- 
ray sources in at least one optical band. Of them, 
322 are detected in the B band, 332 in V, 334 in 
R and 317 in the I band. 

We carried out a quantitative analysis of 
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214.4 214.2 214 

RA (J2000) [deg] 



Fig. 2. — X-ray sources (diamonds) and their opti- 
cal counterparts (filled diamonds) from the cross- 
match of X-ray and optical catalogs 



MAX SEPARATION DISTANCE [arcsec] 

Fig. 3. — Number of optical counterparts as a 
function of the maximum search radius used for 
cross-matching. The number of closest sources 
is represented with stars. Crosses represent the 
total number of sources found within the search 
radius. The dashed line marks the adopted sepa- 
ration. The same value w as obtained applying th e 
statistical methodology of Ide Ruiter et al.l (|1977h 



the completeness and reliabilitjH of our cross- 
m atching procedu r e by a pplying the methodology 
of lde Ruiter et al.l (|l977 ). The separation between 
sources is defined by the dimensionless variable r 



1 Fraction of all real associations between the X— ray and op- 
tical sources that are indeed classified as identifications. 

2 Fraction of claimed identifications expected to be true 
counterparts. 



,Aa' _ 
(-JT + -JT 



af 



1/2 



(2) 



where Aq and A5 are the position differences be- 
tween the X-ray source and its optical counter- 
part in RA & DEC, and a 2 a and a\ are the over- 
all astrometric errors computed as a = (a\ + 
er^) 1 / 2 . In this work, we assume ax =0.7arcsec 
and a pt =0.3 arcsec. ax was calculated by square 
adding a systematic uncertainity of 0.5 and a cen- 
troiding error of 0.5 arcsec. 
The likehood ratio is defined as: 



LR(r] 



dp(r\id) 
dp(r\c) 



-(2X-1) 



2A 



(3) 



where: 



• dp(r|id) is the a priori probability that an 
X-ray source and its optical counterpart 
have a separation between r, and r + dr due 
to astrometrical errors 

• dp(r|c) is the probability that a confusing 
background optical object is found in the 
range r, and r + dr from the X-ray source 
position. 

• A measures the number of confusing ob- 
jects within the error circle and depends on 
the optical objects density (p pt) as: A = 
■na a a$p opt . 

An optical source can be considered as a true 
counterpart of an X-ray source if its LR is higher 
than a certain threshold value L. L is obtained 
maximising the sum of completeness and reliabil- 
ity that are, in turn, mathematically defined in 
terms of the probabilities defined above. In our 
case, L = 0.56 and its corresponding variable is 
r = 2.6 which gives us a separation distance of 2.0 
arcsec, a value consistent with the adopted one. 
The completeness of our detection is 98.3% (i.e. 
some 1 object lost due to the use of a particular 
cutoff value L), while reliability is 87.3% (~ 39 
potentially false matchings). 
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3. Analysis 

3.1. X ray properties 

We have computed the cumulative number 
count distributions (logN-logS) per deg 2 in the 
soft band and compared it with other surveys 
in order check the reliability of our source de- 
tections. We have chosen the soft band to be 
able to co mpare our dist r ibutio n with that com- 
puted by iNandra et al. (2005) and those pro- 
vided by other deep X-ray surveys, Chand ra Deep 
Field South (CDFS, iGiacconi et ahll2002f). Chan- 
dra D eep Field North ( CDFN, I Alexander et al 



2003), ELAIS-N1 and ELAIS-N2 (Manners et al 



2003). All distributions are represented in Fig- 
ure 131 For our sample two curves have been 
computed, one for the complete sample of X— 
ray detections and the other only for the X- 
ray detections that have optical counterparts. 
The median flux levels are 1.81 x 10 -15 and 
1.97 x 10~ 15 ergcm~ 2 s _1 , respectively. For the 
objects brighter than ~ 10~ 14 ergs -1 cm~ 2 both 
distributions coincide, while for objects fainter 
than ~ 10~ 14 ergs -1 cm -2 the density of X-ray 
sources with optical counterparts starts to de- 
crease, since the optical sample is incomplete at 
the depth of our X-ray data. In general, at our 
flux limit the number counts of our complete X— 
ray sample are in good agreement with Nandra 
et al. (2005) and other surveys. We have de- 
tected a higher density of sources with flux below 
~ 5 x 10~ 15 erg cm -2 s _1 , most likely due to the 
use of different detection algorithms and more 
sensitive detection threshold. As we already men- 
tioned, our threshold is 2 x 10 -7 , two times higher 
than the one chosen by Nandra et al. (2005) (see 
12.1.21 for more information). 

Three independent energy bands, soft, hard2 
and vhard (see the Section al. ll for ranges), have 
been used to obtain the X-ray color-color diagram 
shown in Figure [5l Hardness ratios have been 
calculated as described in Section 12.1.21 X-ray 
sources are superimposed on a model grid com- 
puted using power law Spectral Energy Distribu- 
tions (SEDs) with different photon indexes, T, af- 
fected by a galactic absorption column density, 
N#, of 1.3 x 10 20 cm~ 2 and different values of in- 
trinsic absorption; the Chandra ACIS-I response 
matrices have been also included in the computa- 



3.5 r 
3.0 - 
5? 2.5 - 

5T : 

A 2.0 - 
-z. 

cn 
J2 

1.5 - 
1.0 - 




GWS 

GWS (optical counterparts 
GWS, Nandra et al. (2005) 
CDFN, Alexander et al. (2003) 
CDFS, Giacconi et al. (2002) 
ELAIS-N1, Manners et al. (2003) 
ELAIS-N2, Manners et al. (2003) 

-16 -15 -14 

logS (0.5-2 keV) erg cm-2 s-1 



Fig. 4. — Cumulative logN-logS functions for the 
Groth field in the soft band, for all detected X-ray 
sources (thick solid line) and for X-ray sources 
with optical counterpart (thin solid line). Five 
distributions, with different exposure times and 
effective areas, have also been repr e sente d as com- 
parison: GWS by INandra et all (|2005l ) (dotted 
line), CDFN (thick dash-dot-dash line), CDFS 
(thin dash-dot-dash line), ELAIS-N1 (dash line) 
and ELAIS-N2 (dash-dot-dot-dash line) 



tion. To obta in the fluxes we have used XSPEC 
( Arnaud 1996), varying the photon index between 
and 3 with a step size of 0.5, and the intrinsic 
absorption column density between 20 and 24 (in 
logarithmic scale) in steps of 0.5. The majority 
of our sources are located within the model grid 
region, suggesting that their X-ray SED is com- 
patible with the assumed model. A specially high 
concentration of objects can be seen in the region 
of Nh between 10 20 and 5 x 10 21 cm -2 and with 
r~1.8. A group of sources are located outside 
the model grid, showing a soft excess, suggesting 
that they can not be properly described by a single 
power law. For the fitting of these objects proba- 
bly more components should be added. 



3.2. Optical structural parameters 

We obtained a set of structural parameters of 
the optical counterparts of our X-ray sample, in 
order to describe quantitatively their structure 
and morphology. In this section we describe the 
procedure used to derive the parameters and to 
produce a morphological classification. 

We used SExtractor v2.2.2 ( Bertin fc Arnoutsl 
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Fig. 5. — X-ray color-color diagram based on 
2-4.5/0.5-2keV (HR' X ) vs. 4-7/2-4.5keV (HR/,) 
hardness ratios. X-ray sources (crosses), from the 
complete X-ray catalog, are represented and su- 
perposed on an absorbed power law model grid. 
Photon index is varying between and 3 in steps of 
0.5 (dashed lines from up to down) and a log(N#) 
between 20 and 24 in steps of 0.5 (solid lines from 
left to right). Median error bars are represented 
in the upper left corner of the diagram. 



19961 ) to detect and extract the sources. Details 
of the source ext raction process are provided in 
Cepa et al.l (|2008h . The version used (v2.2.2) was 
modified by B. Holwerda to measure the Abra- 



ham concentration index ( Abraham et al.l ll994) 
This is defined as the ratio between the integrated 
flux within certain radius defined by the normal- 
ized radius a and the total flux. We used the val- 
ues obtained at a = 0.3 in our ana lysis in order to 



compa re our results with those of I Abraham et al 
(1996). Another relevant output parameters for 



our study are the CLASSJ3TAR, which permits 
to roughly classify objects as compact or extended 
(for more information see Section [3.2.ip . and the 
background that will be used later. 

In addittion, we used the Galaxy IMage 



2D software pack age (GIM2D, ISimardl [l998; 



Simard et al.ll2002h versions 2.2.1 and 3.1, to per- 



form the galaxy structural parameter decompo- 
sition (see Figure [6]). This software considers 
two components: an spheroid, represented by 
a Sersic profile^, and a disk represented by an 



exponential profi le. Using t he standard steps 
in DAOPHOT (jStetsonl [l987h we obtained the 
PSF images required for GIM2D execution. The 
initial parameters (including background estima- 
tion) and setup for GIM2D were determined using 
the images and catalogs p roduced by SExtrac- 
tor ( Bertin k. Arnoutslll996l ). We obtained from 
the output provided by GIM2D the following pa- 
rameter^: bulge-to-total ratio (B/T), residual 
parameter (R) (S chade et al. 1995, 1996), a sym- 
metry (A) index (jAbraham et al.lll994 [l996h . the 
Sersic index (n) and the total luminosity (L to t)- 
The residual parameter and the asymmetry in- 
dex quantify the galaxy irregularity, by measuring 
the deviation from the simple, symmetric model 
(e.g. spheroid + disk), using as input the residual 
and original image, respectively. GIM2D also cal- 
culates for certain output parameters their lower 
and upper limits at 99% confidence level. In order 
to estimate the accuracy of the best fit parameters 
we have looked at the deviation, computed as the 
maximum between the best value and the upper 
or lower limits (see Figure [7]). In this respect, the 
B/T parameter and Sersic index have revealed in- 
accurate when applied to faint and small objects 
(< 300 pixel, while the seeing disk (FWHM) area 
is about 30 pixels) as shown in Figure[7l Therefore 
we have decided to do not use these parameters 
for the morphological classification. 

In practice, we observed that the code either 
fails or produces unreliable output parameters in 
a number of cases: (1) when the source R magni- 
tude is less than 18 or larger than 24, (2) when the 
isophotal area is very small (usually less then 90 
pixels), (3) when the central object is surrounded 
by close neighbouring sources or nearby bright 
companions, and (4) when the source is close to 



2 While any rigourous analysis of the structural components 
of an AGN host should include an additional component, 



namely a pointlike nucleus, the small size of most of our 
objects compared with the seeing disk prevents from imple- 
menting such a detailed approach. Rather, we have simply 
assumed that the sum of the nuclear and bulge components 
can be roughly represented by a single Sersic profile. 
3 According to the available documentation (Simar d Hl99St 
ISimard et"aLll2002t) , GIM2D calculates the Abraham con- 
centra tion at four normalized rad ii c* = 0.1, 0.2, 0.3 and 
0.4. JAbraham et al.l 11994 1199611 . Comparing the re- 
sults with the ones obtained with SExtractor, a systematic 
shift of ~0.2 is found towards larger values in GIM2D fig- 
ures. After performing several tests we concluded that the 
GIM2D code performs its computations using normalized 
area rather than normalized radius against what is stated 
in the documentation. 
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the frame boundaries. 




Fig. 6. — Four sample images of the GIM2D 
decomposition procedure. First column shows 
subimages (obtained from the scientific image) of 
objects to be fitted; the second column repre- 
sents the mask subimages, third column shows 
the subimages obtained from the modelling and 
the fourth column shows the residual images, ob- 
tained from the original image after subtracting 
the model 
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Fig. 7. — Evolution of the median error of the 
bulgc-to-total ratio B/T (crosses), Sersic index n 
(stars) and total flux (filled circles) l^tot as func- 
tions of the isophotal area. Deviations are evalu- 
ated as explained in Section [ 



3.2.1. Morphological classification 

The selection of compact objects (QSOs, faint 
host AGNs or active stars) has been accomplished 
by using the SExtractor CLASS_STAR output pa- 
rameter. The computed value depends on the 



seeing, the peak intensity of the source and its 
isophotal area. The parameter ranges between 
(extended object like a galaxy) and 1 (pointlikc 
object). A source has been deemed as compact if 
its CLASS_STAR parameter is larger than 0.9. 

To determine the morphology of extended ob- 
jects (i.e. not selected as compact according to 
the criterion above), we analyzed four parame- 
ters: Abraham concentration index (C), Abraham 
asymmetry index (A), B/T ratio and residual pa- 
rameter (R). We carried out the analysis in three 
different bands (B, V, and R) in order to determine 
the most reliable parameters, and the best band 
for morphological classification. We excluded I 
band here due to fringing. 

We first classified each object using the set 
of parameters mentioned above, separately for 
each parameter in each band, considering the re- 
sults from previous works as limits between dif- 
fer ent morphologic al ty pes. We used the results 
of llm et all (|2002h and [SimardJ (|l998l ) to obtain 
preliminary morpho logy with the B/T parameter, 
Schade et al.1 (Il995h resul t s for t he residual param- 



eter and Abraham et al. (1996) for the classifica- 
tion with the concentration and asymmetry index. 
In practice, we found that, in many cases, differ- 
ent parameter sets lead to a different classification 
of a given object even in the same band. In or- 
der to have an independent evaluation criterion, 
we also carried out a visual classification of the 
objects. Despite its subjectiveness, this procedure 
can yield highest-quality results when applied to 
bright and extended objects, thus being a conve- 
nient method to test the reliability of the struc- 
tural parameters that will be in turn applied to 
the analysis of fainter objects. The visual classifi- 
cation has been accomplished by inspection of the 
radial profile, 2D surface plot and isophotal coun- 
tours diagram of each source (see Figure [8]) . For 
the dimmest and smallest objects of our sample 
(galaxies with the magnitude R > 24 and isopho- 
tal area < 100 pixels) the profiles and isophotal 
contours are too noisy to get an acceptable result 
by means of a simple visual inspection. 

From this comparison we conclude that for our 
sample, characterised by faint objects (typically R 
> 22) with small isophotal area (< 300 pixel, while 
the seeing FWHM disk area is about 30 pixels) the 
most reliable parameter is the concentration index 
C, combined either with the assymmetry index or 
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the residual parameter. We therefore based our 
morphological classification in the concentration 
index combined with both asymmetry index and 
residual parameter. The R band has been used, 
since, on the one hand, the largest part of optical 
counterparts have been found in this band, and on 
the other, the achieved S/N is generally larger. B, 
V, and I bands have been used to perform the mor- 
phological classification of objects not detected in 
the R band, or detected in this band but with very 
low S/N ratio. 
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Fig. 8. — Sample of different Hubble type galax- 
ies and their contour diagrams; Up: Left: Group 
I (E/SO), Center: group II (Early type spiral), 
Right: group III (Spiral); Down: Left, group III 
(Spiral), Center: group IV (Irr) , Right: compact 
object (C>0.7 or CLASS.STAR > 0.9). All ob- 
jects are located in the center of the images. 

Based on the classification using visual inspec- 
tion of the brightest sample objec ts, and upon 
comp arison with previous works (e.g. I Abraham et~al 
19961 see references above) in the present paper 
we created five categories or groups associating 
ranges of values of C with Hubble types: 



0: C > 0.7 or CLASS_STAR > 0.9 Compact obj. 
I: 0.45<C<0.7->E, E/SO and SO 
II: 0.3 < C < 0.45 -» SO/SOa-Sa 
III: 0.15<C<0.3^Sab-Scd 
IV: C<0.15->Sdm-Irr 

The first group is populated by elliptical and 
lenticular galaxies, whose profiles do not show any 
trace of disk or spiral arms (see Figure [8]), the sec- 



ond group includes very early spiral type galax- 
ies which show signs of spiral arms. The third 
group comprises spiral galaxies while the fourth 
includes late type spirals and irregular galaxies, 
that show almost no regular form in their profiles. 
All galaxies, even those in the first group (ellip- 
tical and lenticular) have a concentration index 
C < 0.7, while all objects with C > 0.7 have been 
classified by SExtractor as compact ones, with the 
CLASS.STAR parameter > 0.9. 

In Figure [9] we depict the Abraham asymmetry 
vs. concentration index diagram segregated ac- 
cording to the morphology determined using the 
visual classification for the brightest objects, and 
the structural parameters (in most of the cases C 
combined with A) for the fainter o nes. Comparing 



our F igure [9] with that shown in lAbraham et al 
( 19961 ). it can be seen that in our case there is 



a small shift of the concentration index towards 
smaller values. This can be explained consider- 
ing the influence of the atmospheric seeing con- 
ditions. SExtractor does not apply any correction 
on the computed parameters for this effect, and as 
the seeing induced PSF increases, both concentra- 
tion and asymmetry indexes tend to decrease, as 
confirmed by our simulations carried out by con- 
volving simple galaxy mod els with seeing PSFs o f 



(199G) 



different widths. Since the lAbraham et al. 
diagram has been obtained from HST data, these 
differences are expected. 

The morphological analysis described above has 
been applied to the complete sample of 340 X-ray 
objects with optical counterpart (see Figure [TO]). 
Of those, 333 were classified using the R band, 4 
sources using the V band and finally 3 objects by 
means of B band data. The largest group (143 
sources) is composed by compact objects (QSOs, 
faint host AGNs, stars) followed by group III (73 
galaxies). 57 objects are found either in group II 
or III, having outstanding bulge component a signs 
of a formed disk and spiral arms. Only 33 galax- 
ies are found to belong to the IV group. Finally, 
we detected 7 closed pairs, whose possible phys- 
ical relation will be confirmed when photometric 
redshifts are computed (paper III, in preparation). 
Due to the high noise level, or to the presence of 
other objects close to the observed ones, the clas- 
sification uncertainty for some of the objects is 
quite high, and it was not possible to determine 
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Fig. 9. — Morphological classification diagram 
combining the Abraham asymmetry index (A) 
vs. Abraham concentration index (C). Diamonds 
present E and SO galaxies, triangles SO/SOa-Sa 
galaxies, squares young spiral galaxies (Sab-Sbcd) 
and crosses late type galaxies (Sdm-Irr). Com- 
pact objects are represented by star symbols. Dot- 
ted lines separate the locii of elliptical, spiral and 
peculiar galaxie s (from right to left) obtained by 



Abraham et alj (1996). There is a small shift of 
our concentration and asymmetry index towards 
smaller values due to the influence of the atmo- 
spheric seeing conditions. 



the morphological group: 3 galaxies could be ei- 
ther in group II or III and other 3 either in group 
III or IV. For 21 sources it was not possible to 
determine the morphology since they are at the 
detection limit. 

4. Relationship between X-ray and opti- 
cal structural properties 

4.1. Nuclear classification 

The simplest way to classify the nuclear type 
is based on the X-ray-to-optical ratio (X/O). 
The typical value of X/O flux ratios for X- 
ray selected AGN (both type 1 or BLAGN and 
type 2 or N LAGN) is in the range between 0.1 
and 10 (e.g. iFiore et all 120031) . For optically se- 
lected type 1 AGNs th e typical value is X /O ~ 1 
(lAlexander et alJ l200ll IFiore et all l2003h. At 



20011 : IFiore et all l2003h . 
high X/O flux ratios (well above 10) we can 
find BLAGN and NLAGN as well as high-z high- 
luminosity obscured AGN (type 2 QSOs), high-z 
clusters of galaxies and extreme BL Lac objects. 



Fig. 10. — The histogram represents the final mor- 
phological classification of the GWS X-ray emit- 
ters with optical counterparts. Different values of 
the Abraham concentration index correspond to 
different groups I-IV (see text). A group of possi- 
ble interacting objects has been detected; a group 
of sources (marked with '?') remained unidenti- 
fied due to their very low S/N ratio. 

Finally, the X/O < 0.1 region is typically popu- 
lated by coronal emitting stars, normal galaxies 
(both early type and star-forming) and nearby 
heavily absorbed (Compton thick) AGN. 

We have applied a simple criterion for per- 
forming a coarse nuclear-type classification of our 
sample objects, based on diagnostic diagrams re- 
lating X-ray-to-optical ratios (X/O) to hardness 
ratios ( HR). This criterion is b ased on the re- 
sults of iDella Ceca et al.1 l|2004| ) on the XMM- 
Newton Bright Serendipitous Survey; they found 
that, when plotting the BR(2-A.5keV/0.5-2keV) 
(equivalent to our HR^) against X/O, computed 
as the ratio of the observed X-ray flux in the 0.5- 
4.5keV (soft+hard2) energy range and the optical 
R-band flux, most (85%) of the BLAGN identi- 
fied by means of optical spectroscopy are tightly 
packed in a small rectangular region of the dia- 
gram, while NLAGN (type 2) tend to populate a 
wide area of the diagram (towards harder values 
of the hardness ratio). We applied the same ap- 
proach to our sample, combining the HR^ hard- 
ness ratio with the X/O ratio, where the opti- 
cal flux Fr has been derived from the Petrosian 
magnitudes in the R band. The resulting dia- 
gram, that also includes host morphology infor- 
mation, is shown in Figure [TT1 where the dashed- 
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line box indicates the locus of the BLAGN re- 
gion. We set the right edge of this box, HR(2 — 
4.5fceU/0.5 — 2keV) ~-0.35 as a coarse boundary 
separating 'unobscured' AGN (BLAGN) and 'ob- 
scured' AGN (NLAGN). On the other hand, we 
consider the small box depicted in the diagram 
as the "h ighest probability" re gion for finding a 



BLAGN (jDerla Ceca et al l 120041 ). We find that 
a large fraction of our objects, 63% 0, fall inside 
the region of BLAGNs, while 51% of the sample 
is placed within the smal l BLAGN box defined 
bv iDella Ceca et ah (2004). Regarding morphol- 
ogy, we do not find evidence of any relationship 
between nuclear and morphological types, since 
we did not obtain any clear separation of nuclear 
species based on morphological type. Neverthe- 
less, we can see that the majority of objects iden- 
tified as compact (65%) are placed in the BLAGN 
region. It is also interesting to highlight the re- 
lationship between the X/O ratio and the host 
galaxy morphology, with early-type galaxies hav- 
ing generally lower values for X/O ratio and an 
opposite behavior of late-type objects. Finally, 
we find that only about 7% of our sample objects 
are placed in the lower region of the diagram (i.e. 
that typical of stars/normal galaxies/Compton- 
thick AGN). 

We also represented the X/O ratio with respect 
to HRi, yielding a basically identical distribution 
as that described above, and with respect to HR2 
and HR' 2 , observing that in these two cases is 
harder to get a clear separation between differ- 
ent nuclear/morphological types; this behavior is 
expected since hardness ratios computed from the 
hard and vhard bands are less sensitive to absorp- 
tion than those involving the soft band, which is 
one of the main criteria for the separation between 
broad and narrow line AGN. 

4.2. X-ray and optical properties 

We combined the X/O ratio with optical col- 
ors (B-R and B-I) for the different morphologi- 
cal groups and nuclear types, as shown in Fig- 
ure [12] We observe that compact and late-type 
objects tend to show bluer colors, while early- 



4 A11 fractions given here are relative to the total number 
of X— ray objects having optical counterpart and detected 
in both soft and hard2 energy bands such as HR(2 — 
4.5fceVy0.5 — 2keV) can be computed 



type galaxies tend to present redder colors and 
also lower X/O values, as expected. On the other 
hand, there is a mild tendency for BLAGN to show 
bluer colors than those found in NLAGN, as al- 



ready objerrod_2rij£raljaniveree_samples (e.g.|Yee 



bay 

1983 ISanchez-Portal et al.ll2004h . The B-R av 



erage value for BLAGNs is 1.0 and for NLAGNs 
1.4, with a — 0.5 and 0.6 respectively. Using the 
Kolmogorov-Smirnov test, the two distributions 
are statistically different with a significance level 
of 99%. 

When representing the X/O ratio vs. the Abra- 
ham concentration index C we find a strong anti- 
correlation between these two parameters (Figure 
[15]) , obtaining a correlation coefficient of ^-0.5. 
Using Spearman and Kendall (S-K) statistics, we 
obtained a correlation significance level >99%. 
The first order polynomial function has been fit- 
ted with a slope of -2.03, logX/O-intercept of -0.94 
and a standard deviation of ^0.5. This result is 
consistent with that referred above and depicted 
in Figure [Til t na t is the tendency of early-type 
galaxies to have lower than average X/O values 
and the opposite behavior in late-type ones. As 
seen in Figure 1131 all morphological groups, in- 
cluding compact objects, tend to follow the same 
relation. For compact objects we obtained simi- 
lar parameters using S-K statistics: a correlation 
coefficient of ^0.5 and a correlation significance 
level >99%. A least squares fit gives to the anti- 
correlation a slope of -2.4 and a logX/O-intercept 
of -0.95. Comparing different nuclear types, we 
found that both BLAGN and NLAGN follow the 
same distribution (see Figure IT3|> . On the other 
hand, we did not find any clear correlation be- 
tween HR(2 - 4.5A:eU/0.5 - 2keV) and the con- 
centration index (Figure [14]) , suggesting that the 
observed anticorrelation is not related to obscura- 
tion. 



The physical origin of the observed anticorre- 
lation is not well established. First of all, it is 
not easy to disentangle possible bias effects that 
superimposed to variable seeing conditions. It 
has been pointed out the existence of redshift- 
dependent biases in both A and C, due to the 
reduction of apparent size and surface brightness 
with respect to the sky background with increas- 
i ng z, along with poorer s ampling and lower S/N 
(Brinchmann et al. 1998; IConselice et al. 2003); 
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Fig. 11. — Relationship between X/O ratio and HR(2 — 4.5fceV/0.5 — 2keV), for the different morphological 
types. Solid lines separate BLAGN and NLAGN r egions and the are a with X/O < 0.1, where stars, normal 
galaxies and Compton thick AGNs can be found ( Fiore et al.ll2003h . Dashed line box presents the limits 
obtained by iDella Ceca et~aL (2004) where ^85% of their spectroscopically identified BLAGN have been 
found. In our case, inside of this box can be found ~ 51% of all objects from the diagram. 
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Fig. 12. — Log(X/0) flux ratio vs. B-R color for different morphology (left) and nuclear types (right). 



on the other hand, there is a "ban dpass shift" with 



respe ct to the rest wavelength (|Bershady et al 
2000). Even if these effects combine to produce 



some uncertainty in the morphological classifica- 
tion of galaxies, generally in the sense of shift- 
ing objects to apparently later Hubble types (i.e. 
with lower C values), A and C seems to be repro- 
ducible out to z = 3, within a reasonable scatter 
( Conselice et al.l [20031 ). Moreover, there can be a 
truly evolutionary effect in the AGN host galaxies, 
indicating that galaxies at hi gher redshift are i n- 
trinsically less concentrated (jGrogin et al.ll2005l ). 
At this stage we cannot apply any correction 



which depends on the redshift, although this is 
one of the objectives for the future work. How- 
ever, we have performed a preliminary evaluation 
by representing the concentration index C versus 
the spectroscopic redshifts z gathered from the 
DEEP2 catalog for 84 objects from our sample. 
While for non compact objects we obtained an an- 
ticorrelation between these two parameters, C and 
z are found uncorrelated for compact objects, in 
clear contrast with the observed behavior in Fig- 
ure m 

So, although the biases described above might 
contribute partially to it, the observed anticorre- 
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Fig. 13. — Relationship between X/O and the Abraham concentration index C for different morphological 
(left) and nuclear types (right). A clear anticorrelation is observed. 
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Fig. 14. — Relationship between the hardness ratio HR(2 — 4.5keV/0.5 — 2keV), and the concentration index 
C for the different morphological types. 



lation seems to also be a result of a true physical 
effect. In order to attempt to explain it, we have 
first to understand the physical connections of the 
two observables C and X/O. The concentration 
index (C) is tightly linea rly correlated to the c en- 
tral velocity dispersion ( Graham et al. 2001 al ibi), 
which is in turn related to the nuclear b lack hole 
(BH) mass (e.g. iMagorrian et al.l ll998). There- 



fore, the concentration and the nuclear BH mass 
(Mbh) are correlated. The physical connections 
of the X-ray to optical flux ratio are not so ev- 
ident. The X/O ratio measures the X-ray flux 
(in the 0.5-4.5keV band), normalized with the R- 
band optical flux of the whole galaxy (nucleus, 
bulge and disk) . It should be noted that no band- 



pass, or K-correction is performed; therefore if the 
sample spans a large redshift range the interpre- 
tation of the quantity becomes more problematic. 
If the nuclear luminosity is large compared with 
that of the host galaxy, X/O ratio can be thought 
as a measure of the X-ray to optical spectral in- 
dex, aoJE This parameter has been found to 
be strongly anticorrelated with the UV luminos- 
ity L ° , without a significant correlation with 

J 2500i5r! — 1 

redshift (Stc ffen et al. 



20061) . On the other hand 
Bian (2005) found a strong correlation between 
the hard/soft X-ray spectral index, ax and the 



3 This relation is usually defined as aox 
0.3838 log(/ 2fee W/ 250oA ) 
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Eddington ratio in a sample of 41 BLAGN and 
narrow-line Seyfert 1 (NLS1) galaxies observed 
with ASCA . Other wor k s have confirmed thi s 
correlation (jGrupel 12004 IShemmer et~ai1 120061 ). 
The relationship between ax and aox has not 
been thoroughly studied in large samples of AGN, 
but a linear correlation has been found between 
both spectral indexes in a sample comprising 22 
out of 23 quasars in the complete the Palomar- 
Gre en X-ray sample with z<0.4 and Mb <- 



23 ( Shang et al. 20071 ). If this last correlation 



holds for our sample, X/O could be tracing the 
Eddington ratio in the large nuclear luminosity 
limit. In addition, if the host galaxy luminosity is 
large when compared with the nuclear luminosity 
(Lbuige+disk > L nuc i eus ), the X/O ratio can be 
thought as a lower bound of the X-ray-to-bulge 
luminosity ratio, that is in turn a (weak) measure 
of the AGN Eddington ratio L/LEdd oc L/Mbh, 
assuming that the X-ray luminosity represents 
the nuclear luminosity and the bulge luminosity is 
proportional to the bulge mass (and therefore to 
M BH ). 

Hence, we can guess a (loose) correlation be- 
tween X/O and the energy production efficiency 
of the AGN. Under these assumptions, the C 
vs. X/O relation traces the correlation between 
the Eddington ratio and the nuclear BH mass. 
The obtained result could therefore suggest that 
early-type galaxies, having poor matter supply to 
feed the activity, have lower Eddington rates than 
those of late-type, with larger reserves of the gas 
for AGN feeding. 

This suggested appr oach is consisten t with 
the results obtained by iBallo et all (|2006l . l2007h 
in a sample of X-ray selected AGNs at z<l. 
They observed that a large fraction of low- 
luminosity AGNs are fed by black holes with 
a mass Mbh > 3 x 10 6 M©, and Eddington rate 
< 1, while a small ( < 10%) fraction of AGNs 
host small SMBH (M BH < 10 6 M©) and Edding- 
ton rate ~ 1, which corresponds to a more effi- 
cient accretion rate. They stated that their results 
strongly suggest that most of the low-luminosity 
X-ray selected AGNs (at least up to z < 0.8) are 
powered by rather massive black holes, experienc- 
ing a low-level accretion in a gas-poor environ- 
ment. 

Our approach i s also consistent with the re- 
sults of IWu fc Liul (|2004h . who studied the black 



hole masses and Eddington rates of a sample of 
135 double-peaked broad line AGNs, observed in 
two surveys: SDSS and a survey of BLAGN ra- 
dio emitters. They obtained that if the separa- 
tion between the line peaks (which is correlated 
to the SMBH mass) decrease, the Eddington rate 
increases. 

iKawakatu et al. ( 2007 ) found an ant icor relation 
between the mass of a SMBH and the luminosity 
ratio of infrared to active galactic nuclei Edding- 
ton luminosity, L/^/L^d, for a sample of ultralu- 
minous infrared galaxies with type 1 Seyfert nuclei 
(type 1 ULIRGs) and nearby QSOs, which also 
could be consistent with our approach. The an- 
ticorrelation is interpreted as a link between the 
mass of a SMBH and the rate of mass accretion 
onto a SMBH, normalized by the AGN Eddington 
rate, which indicates that the growth of the black 
hole is determinated by the external mass sup- 
ply process, and not the AGN Eddington-limited 
mechanism, changing its mass accretion rate from 
super-Eddington to sub-Eddington. 

5. The clustering of X-Ray emitters and 
extended optical counterparts 

The study of large-scale environmental condi- 
tions of AGN host galaxies are useful not only 
for understanding the differential properties of the 
AGN phenomena, but to trace the formation and 
evolution of such galaxies and the structures where 
they reside. Despite the relevance of X-ray se- 
lected AGNs, their cl ustering properties sti ll re- 



main quite unknown (Basilako s et al. 2005, and 



references therein). In this sense, and taking into 
account the characteristics of the X-ray source 
catalog presented in this work, we have calculated 
the 2-point angular correlation function (2p-ACF) 
for significant subsamples. 

The 2p-ACF, u>(6), gives the excess of probabil- 
ity, with respect to a random homogeneous distri- 
bution, of finding two sources in the solid angles 
dfli and dD,2 separated by an angle 9, and it is 
defined as 



dP=JV 2 [l + u(0)} 



(4) 



where M is the mean number density of cataloged 
sources. To measure the 2p- ACF we have follow 



the same procedure as do ne in|C cpa et al 



adopting the estimator of lLandv &: Szalav 



(2008), 



1993), 



1G 



which can be written in the form 

DD -2DR + RR 



RR 



(5) 



where DD — N ss (9 i ) is the fraction of possible 
source-source pairs counted in z-bins over the an- 
gular range studied and DR = 



(N g -1) 
2N r 



N gr (9i), 



RR 



N S (N S ~1) 



N rr (9i) are the normalized 
random-random 



Nr(Nr-l) 

counts of source-random and 
pairs in these bins, respectively. To make the ran- 
dom catalogs we placed (200 x N source ) random 
points by following the X-ray source angular dis- 
tribution. Howewer, the source density is affected 
by several instrumental biases such as detector 
gaps and sensivity variations due to survey re- 
gions with different exposure times. To account 
for these effects, we made normalized sensitivity 
maps from the effective exposure time and source 
distributions for each band used in this section 
(soft and hard2). The raw maps were modelled 
using a mimimum-curvature algorithm. Next, the 
uniform density random catalogs were convolved 
with the corresponding sensitivity map to obtain 
differential number counts in the surveyed area. 
In the case of the HR(2 - 4.5fceV/0.5 - 2keV) 
catalog, we adopt the sensitivity map used for the 
hard2 source one. The smoothed footprint of each 
sensitivity map can be seen in Figure 1151 

We measure w(9i) in i— logarithmic bins of 
width A (log 9) = 0.2 in the angular scales 0.24 < 
9 < 24.0 arcmin. The lower limit in angular sepa- 
ration is according to the minimum detection scale 
defined in the Section [2. 1.21 The formal error as- 
sociated with the 2p-ACF measured is the Poisso- 
nian estimate ct^ = . The 2p-ACF binned 

measures were fitted using the A UJ (9 1 ~ 1 — C) power 
law. We adopted a fixed canonical slope 7=1.8 
(the so called "comoving clustering model") and 
C is defined as 

c=^m^ (6) 



E N rr (9) 

corresponding to the nu merical estim ation of 
the "integra l const raint As in 

Mivaii et al.l (|2007l ). we use the normalization A 
as the fitting parameter rather than the correla- 
tion length 9 C — A 1 /(~'~ 1 \ since the former gives 
better convergence of the fit. 

From the whole catalog we have selected the 
soft (0.5-2 keV) and hard2 (2-4.5 keV) sources 



soft (0.5-2 keV) 



214.8 214.6 214.4 214.2 214.0 213.8 213.6 213.4 

Right Ascension [degrees] 



hord2 (2-4.5 keV) 




214.8 214 



214.4 214.2 214.0 213 
Right Ascension [degrees] 



213.6 213.4 



Fig. 15. — Smoothed sensitivity maps in the soft 
(0.5-2 keV) and hard2 (2-4.5 keV) bands, used 
to convolve the uniform density random catalogs. 
The sensitivity dimension is normalized to the 
maximum exposure time in each band. Contours 
are linearly binned and lighter regions correspond 
to higher sensitivities. 



(given their similar depths in the survey), the X- 
ray objects with HR(2 - 4.5keV/0.5 - 2keV)^0 
and the fraction of X-ray sources with optical 
counterparts. Table [5] summarizes the basic re- 
sults of the correlation analysis for these sub- 
samples. In all cases, we have obtained signif- 
icant positive clustering signal. Regarding the 
soft band sample, the correlation length 9 C de- 
rived from is cons istent with the obtained by 
Basilakos etafl (|2005h (XMM-Newton LSS) and 
it is ~2a la r ger t ha n the values o btaine d by 
Gandhi et all (120061) . ICarrera et all (|2007t ) and 
Ueda et al. ["~ (|2008f) . In the specific case of the 



hard2 band, we only found in the li teratu re the 2p- 
ACF estimations of iMivaii et all (|2007l) . valid for 
X-ray sources detected in the XMM-Newton ob- 
servations of the COSMOS field and correspond- 
ing to their MED band: our fitted correlation 
length is ^4.3 times larger, although our sample 
is a factor 2 deeper. 

In Figure [TH] are represented our 2p- ACF mea- 
surements with the best-fit A u 9~ o s power-law 
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Tabic 2: Clustering properties of the X-ray selected sources. N source is the number of data points used in 
the 2p-ACF estimation and A u is the best- fit value in the separation range given by # m ; n — max , assuming 
a power-law u>(9) — 0( 1_ 7)=-°- 8 , The goodness-of-fit statistics (x 2 ) and the fitted 2p-ACF at 1 arcmin 
[w(l')] are given in the final columns. The errors represent la uncertainties. 



Subsample 


^source 


^min ^max 


A, (1-7 = -0.8) 


x 2 


c(l') 






arcmin 


xl0~ 3 degrees 






soft (0.5-2 keV) 


465 


0.24-24.0 


6.99 ±0.68 


0.030 


0.15 ±0.02 


hard2 (2-4.5 keV) 


266 


0.24-24.0 


9.19 ± 1.28 


0.107 


0.20 ±0.03 


HR(2-4.5fceV/0.5-2fceV) 


190 


0.24-24.0 


9.68 ± 1.30 


0.109 


0.21 ±0.03 


Optical counterparts (I<23) 


104 


0.60-24.0 


19.99 ±3.32 


0.096 


0.42 ±0.07 



in soft and hard2 bands (panels a and b, respec- 
tively). It is noticeable that the clustering ampli- 
tude obtained in both bands looks similar. Thus, 
in the case of the Groth field, we would not expect 
redshift or spatial clustering differential effects in 
the soft and hard2 band samples. Interestingly, 
an absence of correlation in diffe rently defined 
hard bands has been repo r ted b y iGandhi et al 



(2006) and iPuccetti et a 
Yang et all {2003J) and" 



(120061): never t heless , 



Basilakos et all (2004) 



not only yield a positive correlation amplitude, 
but the hard band signal is several times that 
of the soft band. Other results with similar 
correlation le n gths i n both bands ar e giy en by 



Yang et all |2006) 



Carrera et alj (120071) . iMivaii et all (|2007f ) and 



In the panel c of FigurefTBlis represented the 2p- 
ACF estimations fo a HR(2 - 4.5keV/0.5 - 2keV) 
selected sample. The clustering amplitude is 
slightly larger than the soft or hard2 band ones. 
Moreover, the whole HR(2 - 4.5feeF/0.5 - 2keV) 
subsample was decomposed using the bound- 
ary referred in Sect. 14.11 to separate type 1 (or 
BLAGN) from type 2 (or NLAGN) populations 
and a positive, very similar clustering signal from 
both subsamples emerged from the 2p-ACF es- 
timations, despite their goodness-of-fit statistics 
are worse than the tabulated values in a factor 
~ 1 because of the small size of the subsamples. 

Finally, we measured the 2p-ACF of the op- 
tical counterparts of X-ray sources in order to 
compare the clustering behavior with those that 
come from the analysis of galaxy populations in 
the same field. For a reliable comparison is nec- 
essary to diskard the compact objects (i.e. with 
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Fig. 16. — 2p-ACF estimates of the X-ray sources 
detected in the soft band (panel a), hard2 band 
(panel b) and HR as defined in eq. 1 (panel c). 
The continuous line represents the A w 0( 1_ 7)=-O- 8 
power law fit to the data. 

CLASS_STAR>0.9) from the optical counterparts 
catalog, as well as all objects fainter than a lim- 
iting magnitude. We have selected the peak of 
the apparent brightness distribution in the I-band 
as such limiting magnitude. Thus, X-ray coun- 
terparts with I>23 were ruled out, remaining 
104 extended objects, which can be regarded as 
AGN host candidates of the corresponding X-ray 
sources. 
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In spite of this skinny sample, it is possible 
to obtain some interesting lines of thought about 
the X-ray selected AGN clustering. It is sug- 
gested that the typical environments of AGNs 
are no different to those of inactive galaxies 
in general, at red s hift ranges above z~ 0.4 (e.g. 



Grorin et all 12003 



thcless, ICoil et al 



.Waskett et alJl2005h. Never - 
(|2004T ) and iGilli et alJ (|2005h 



have argued that the AGN are preferentially 
hosted by early-type gal axies at a typic a l red - 
shift of ~ 1 . In this way, iBasilakos et al. (2004) 
found that clustering length values of hard X-ray 
objects are comparable with those of Extremely 
Red Objects ( EROs) and luminou s radio sources. 
More recently, iGeorgakakis et al" ( 2007h demon- 
strated that X-ray-selected AGNs at z~ 1, in the 
AEGIS field, avoids underdense regions at 99.89% 
of significance. 

In Figure [17] we have depicted the 2p-ACF 
of the optical counterparts selected as described 
above. The 2p-ACFs of the Groth field galaxy 
samp le with K23 only and with V-I>3 (I<24) 
from ICepa et al.l (|2008l ) are overplotted. At least 
in the range 0.01 < <0.1 degrees, the cluster- 
ing signal produced by the optical counterparts of 
our X-ray sources is l-2<r stronger than the corre- 
sponding to the I< 23 galaxy population. In terms 
of fitted correlation length, the 9 C (X-ray) (1-7=- 
0.8) is ~ 30 times larger than that of the Groth 
field galaxies, but only ~ 2.3 times larger than the 
maximum correlation length of the V-I>3 galax- 
ies in the same field. In the other hand, using 
the optical data from the counterparts catalog, we 
found that the median V-I observed color is 1.86 
and most of 65% of the optical counterparts have 
an observed color V-I>2. This means that a sig- 
nificant population of our optical counterparts is a 
composition of red and very red galaxies, whereas 
its clustering signal tends to reproduce that asso- 
ciated with such galaxies, which driv e the ~ 10x 
clustering excess found by ICepa et al.l (|200S ) . At 
once, more than a 55% of the optical counterparts 
belong to the morphological categories I and II, as 
defined in the Section f3. 2. II 

6. Summary and Conclusions 

We processed and analyzed public, 200 ksec 
Chandra/ ACIS observations of three fields com- 
prising the original GWS gathered from the Chan- 




log({5) [degrees] 

Fig. 17. — 2p-ACF estimation for the opti- 
cal extended (CLASS_STAR<0.9) counterparts 
with I magnitude brighter than 23 (solid cir- 
cles). The continuous thick line represents the 
A,, 0(i-t)=-o.8 power law fit t0 the data ^ The 

dashed line represents the 2p-ACF estimation for 
galaxies with I<23 in the Groth field, whereas the 
continous thin line it does for the V-I>3 color se- 
lected galaxies w ith I < 24, both extracted from 
Cepa et al.1 ([20081 ). The shaded regions represent 
fitting errors of la. 

dra Data Archive, combined with our broadband 
survey carried out with the 4.2m WHT at La 
Palma. Our X-ray catalog contains 639 unique 
X-ray emitters, with a limit flux at 3tr level of 
4.8 x 10 -16 erg cm -2 s _1 and a median flux of 1.81 
x 10~ 15 ergcm~ 2 s _1 in the full (0.5-7 keV) band. 
Cross-matching the X-ray catalog with our opti- 
cal broadband catalog, we have found 340 X-ray 
emitters with optical counterparts. The complete- 
ness and the reliability of cross-matching proce- 
dure are 97% and 88.9%, respectively. At our 
flux limit, the number counts of our complete X- 
ray sample are in good agreement with Nandra 
et al. (2005) and other comparison surveys men- 
tioned in Section 13.11 X-ray SEDs of the mayor- 
ity of our sources are compatible with the power 
law model, affected with the galactic absorption 
column density and the intrinsic absorption. We 
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have performed a morphological classification ob- 
taining and analyzing different structural param- 
eters, including the bulge-to-total ratio, residual 
parameter, Sersic index, asymmetry and concen- 
tration indexes. We also performed the classifica- 
tion based on visual inspection. To obtain the nu- 
clear classification, we applied a criterium based 
on the correlation between X/O flux ratio and 
HR1' hardness ratio. 63% of our sources, with 
HRl'^0, have been classified as BLAGNs. The 
electronic table summarizes the X/O flux ratios, 
HR1' values and morphological classification of all 
340 X-ray emitters with optical counterparts. We 
have investigated the correlations between X-ray 
and broadband optical structural parameters that 
provide with information about the host galaxy 
morphology and populations, finding that: 

1. For a sample where most of the objects are 
faint (R > 22) and with small isophotal area (< 300 
pixels), bulge-to-total flux ratio has revealed inac- 
curate to perform the morphological classification. 
Concentration index combined with the asymme- 
try index have been found to provide the best 
separation criterion for such a sample, different 
morphological types have been observed to dis- 
tribute according to the following ranges of con- 
centration index: group (compact objects) — > 
C > 0.7 (or CLASS.STAR > 0.9), group I (E, E/S0 
and SO) -> 0.45<C<0.7, group II (SO/SOa-Sa) 
-> 0.3<C<0.45 (0.47), group III (Sab-Scd) -> 
0.15 < C < 0.3 and group IV (Sdm-Irr) C < 0.15 

2. No clear separation has been found between 
the morphology and the nuclear type, but most of 
our objects classified as compact are placed in the 
region of BLAGNs. 

3. We can confirm the tendency of compact and 
late type objects to show bluer colors, while early 
type galaxies tend to present redder colors. On the 
other hand, we find a mild tendency of BLAGNs to 
show bluer colors when compared with NLAGNs. 

4. We obtained an anticorrelation between the 
X-to-optical flux ratio and the concentration in- 
dex (which parametrizes the morphology), with 
significance level higher than 99%, showing that 
the early-type galaxies tend to have lower and the 
late- type galaxies higher than average X/O values. 
Objects classified as compact show this anticor- 
relation as well, with the same significance level. 
The anticorrelation might suggest that early-type 
galaxies, having poor matter supply to feed the ac- 



tivity, have lower accretion rates than those of late 
type, with larger reserves of gas for AGN feeding. 

Finally, we presented the first results on the 
two-point angular correlation function of the X— 
ray selected sources in the Groth field from Chan- 
dra /ACIS observations, segregated in four sub- 
samples: Soft, hard2, HR( 2 _4.5fc e y/o.5-2fceV) an d 

optical extended counterparts. We obtain the fol- 
lowing results and conclusions: 

1. All subsamples reveal a significant positive 
clustering signal in the ~0.5 to 6 arcmin separa- 
tion regime. Canonical power-law (7=1.8) fits to 
the angular correlation functions have been made, 
obtaining correlation lengths between ~ 7 and 27 
arcsec, depending on the subsamplc. 

2. We have neither found substantial differ- 
ences between the clustering of soft and hard2 
selected X-ray sources, nor when comparing the 
populations of type 1 and type 2 AGN (i.e. 'unob- 
scured' and 'obscured' AGN, respectively), sepa- 
rated using the HR(2-4.5/ceF/0.5-2fceF) w - 0.35 
boundary. 

3. The clustering analysis of the optical coun- 
terparts of our X-ray emitters provides a correla- 
tion length significantly larger than the one found 
for the whole galaxy population with the same lim- 
iting magnitude, and same field, but similar to 
that of strongly clustered populations like red and 
very red (bulge-dominated) galaxies. A clustering 
signal > lcr was observed even at ~ 30 arc- 
sec, which corresponds to 160 hr 1 kpc, assuming 
a mean redshift of 0. 85 for Chandra extragalactic 
sources (|Hasanll2007h . This comoving distance is, 
for instance, comparable to the separation of two 
mergers in the t«2 Gyr stage of the ir large-scale 
evolut ion, according to the model of iMaver et al 



(|2007l ) on rapid formation of SMBH. Hence, our 
results suggest that the environment plays an im- 
portant (but not unique) role as possible triggering 
mechanism of AGN phenomena. This fact can be 
perfectly compatible with the interpretation given 
to the anticorrelation between X/O ratio and the 
concentration index. The former may be associ- 
ated to the genesis of the AGNs, while the anti- 
correlation could be that to their maintenance and 
evolution. 
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ABSTRACT 

The Groth field is one of the sky regions that will be targeted by the OTELO (OSIRIS Tunable 
Filter Emission Line Object) survey in the optical 820 nm and 920 nm atmospheric windows. This 
field has been observed by AEGIS (All-wavelength Extended Groth strip International Survey) 
covering the full spectral range, from X-rays to radio waves. Chandra X-ray data with total 
exposure time of 200ksec are analyzed and combined with optical broadband data of the Groth 
field in order to study a set of structural parameters of the X-ray emitters and its relation with X 
ray properties. We processed the raw, public X-ray data using the Chandra Interactive Analysis 
of Observations and determined and analyzed different structural parameters in order to produce 
a morphological classification of X-ray sources. Finally, we analyzed the angular clustering of 
these sources using 2-point correlation functions. We present a catalog of 340 X-ray emitters 
with optical counterpart. We obtained the number counts and compared them with AEGIS data. 
Objects have been classified by nuclear type using a diagnostic diagram relating X-ray-to-optical 
ratio (X/O) to hardness ratio (HR). Also, we combined structural parameters with other X-ray 
and optical properties, and found for the first time an anticorrelation between the X/O ratio 
and the Abraham concentration index which might suggest that early type galaxies have lower 
Eddington rates than those of late type galaxies. A significant positive angular clustering was 
obtained from a preliminary analysis of 4 subsamples of the X-ray sources catalog. The clustering 
signal of the opticaly extended counterparts is similar to that of strongly clustered populations 
like red and very red galaxies, suggesting that the environment plays an important role in AGN 
phenomena. 

Subject headings: Galaxies: active — X-rays: galaxies 



1. Introduction 

The OTELO project (Cepa et al. 2005, 2008) 
is a flux-limited survey of emission-line objects 
in large and perfectly defined volumes of the 
Universe. The OTELO survey, more than one 
magnitude deeper than other emission line sur- 
veys, includes a wide variety of astronomical 
sources: star-forming galaxies, starburst galax- 
ies, emission-line ellipticals, AGNs, QSOs, Lya 
emitters, peculiar stars, etc. With these data, a 
number of scientific problems will be addressed, 
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including: star formation rates (SFR) and metal- 
licity evolution of galaxies; AGN and QSO evolu- 
tion. Deep imaging programs such as several VLT 
surveys for GOODS (Dickinson et al. 2004), COS- 
MOS (Capak et al. 2007) and SXDS (Furusawa 
et al. 2008) fields will serve as preparatory work 
for OTELO. Likewise, we have performed a deep 
BVRI broadband survey in one of the fields of the 
Groth strip selected for the OTELO survey (Cepa 
et al. 2008). 

In the present work, data from this broad band 
survey is combined with other wavelength data 
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for achieving several objectives of the OTELO 
project. Specifically, the matching of these op- 
tical data with X-ray allows succesfully tackling 
the study of the AGN population as shown in Bar- 
cons et al. (2007), Georgakakis et al. (2006), and 
Steffen et al. (2006). Optical imaging surveys al- 
low a morphological classification of all sufficiently 
resolved sources. This is essential to reveal the na- 
ture of AGN host galaxies (e.g. Hasan 2007; Fiore 
et al. 2001), although the morphological classifi- 
cation becomes difficult at high redshift, due to 
the reduced resolution of images. Several works 
(Graham et al. 2001b; Ferrarese & Merrit 2000; 
Gebhardt et al. 2000) have shown that the AGNs 
are directly related with some host galaxy proper- 
ties, particularly with their bulges. Kauffmann et 
al. (2003) have pointed out that the AGNs in the 
local universe are hosted predominantly in bulge- 
dominated galaxies. More recently, Hasan (2007) 
using HST and Chandra data of the GOODS 
South field has found that the most moderate lu- 
minosity AGNs hosts are bulge-dominated in the 
redshift range z~ 0.4-1.3. A tight correlation be- 
tween black hole mass and the bulge dispersion ve- 
locity has been well established (Ferrarese & Mer- 
rit 2000; Gebhardt et al. 2000), as well as with the 
concentration of bulges (Graham et al. 2001b). 

On the other hand, deep X-ray extragalactic 
surveys are extremely useful tools to investigate 
the time evolution of the population of active 
galactic nuclei (AGN) and to shed light on their 
triggering mechanisms. The exceptional effective- 
ness at finding AGN arises largely because X-ray 
selection has reduced absorption bias, minimal di- 
lution by host-galaxy starlight and allows efficient 
optical spectroscopic follow-ups of high probabil- 
ity AGN candidates with faint optical counter- 
parts. The Chandra and XMM- Newton observato- 
ries have revolutionized this field of research, gen- 
erating the most sensitive X-ray surveys ever per- 
formed. In fact, deep surveys produced by these 
observatories have increased the resolved fraction 
of the cosmic X-ray background (CXRB) to about 
90% in the 0.5-2 keV range. A large part of the de- 
tected X-ray sources are AGN (< 70%), although 
X-ray emitters include a mix of different types 
of objects such as galaxy clusters and stars. On 
the other hand, these powerful observatories make 
it possible to access the hard (2-10 keV) range, 
allowing to directly probe the AGN activity not 



contaminated by star formation processes. More- 
over, hard X-ray surveys are capable to detect all 
but the most absorbed (Compton thick) sources. 
Therefore, they provide the most complete and un- 
biased samples of AGNs (Mushotzky 2004; Brandt 
& Hasinger 2005). Number counts relations have 
now been determined (e.g. Cappelluti et al. 2007) 
although with some evidence for field-to-field vari- 
ations. Such variations are expected at some level 
owing to " cosmic statistics" associated with large- 
scale structures and clustering features that have 
been detected in the X-ray sky (Barger et al. 2003; 
Yang et al. 2003; Gilli et al. 2003, 2005). 

In this paper we present the analysis of pub- 
lic, deep (200 ksec) Chandra/ ACIS observations 
of three fields comprising the original Groth- 
Westphal strip (GWS), gathered from the Chan- 
dra Data Archive, and combined with optical 
BVRI data from our broadband survey carried 
out with the 4.2m William-Hcrschel Telescope 
(WHT) at La Palma (see Cepa et al. 2008, for 
a detailed description). We have tried to imple- 
ment an innovative approach to the optical study 
of X-ray emitters by investigating correlations 
between broadband optical and X-ray parame- 
ters, specifically the X-ray-to-optical ratio and 
hardness ratios, and optical structural parameters 
that provide information about the host galaxy 
morphology and populations, like asymmetry and 
concentration indexes and optical colors. 

This paper is structured as follows: in Section 
2, we describe the observational data, including 
X-ray data processing and source detection, as 
well as optical broadband data and the selection 
of the sample of X-ray sources with optical coun- 
terparts. Section 3 reviews the X-ray properties 
of the galaxy sample, the procedures for comput- 
ing optical structural parameters, as well as the 
morphological analysis carried out. The reliabil- 
ity of different structural parameters used in mor- 
phological classification is also discussed. Section 
4 presents the nuclear type classification, the re- 
lationship between X-ray and optical structural 
properties, and the comparison with results from 
the literature. The catalog of X-ray emitters with 
optical counterparts and their morphological clas- 
sification is presented in the electronic edition. Fi- 
nally in Section 5 we present the results of the 
clustering analysis (based on the angular corre- 
lation function formalism) applied on 4 selected 
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subsamples of the X-ray catalog. 

2. Observational data 
2.1. X ray data 

2.1.1. Observations and data processing 

Chandra has observed three consecutive fields 
centered at the original HST Groth-Westphal 
strip (GWS) using the ACIS-I instrument. All 
datasets have been gathered from the Chandra 
Data Archive (http://asc.harvard.edu/cda/) 
using the Chaser tool. PI of all the retrieved 
Chandra observations is K. Nandra. Total expo- 
sure time in each field is about 200 ksec. The total 
field size of ACIS-I chips 0, 1, 2 & 3 is 16.9' x 16.9' 
(ACIS-S chips 2 & 3 were also used, but their 
FOVs do not overlap with our optical data and 
therefore were not used). The FAINT and VFAINT 
telemetry modes have been used, with telemetry 
saturation rates of 170 and 67 events/s, respec- 
tively. Due to the lack of bright sources within 
the fields, we assume that there are no pile-up ef- 
fects in any observation. Table 1 summarizes the 
main characteristics of the Chandra observations. 

We processed the data Using the Chandra 
Interactive Analysis of Observations (CIAO) 
(http : / /exc . harvard, edu/ ciao/), v4.0 and Cal- 
ibration Data Base (CALDB) v3.4.0. Standard 
reduction procedures have been applied, as de- 
scribed in the CIAO Science Threads to produce 
new level 2 event files. The output level 2 event 
files have been restricted to the 0.5-8 keV range 
to avoid high background spectral regions. Fur- 
thermore, we analyzed the source-free light curves 
in order to define additional good time intervals 
(GTI) . To this end, we applied the CIAO detection 
program celldetect with a threshold S/N = 3. 
The detected sources have been removed upon 
creation of the time-binned light curves, that have 
been then used for the computation of the GTIs. 
A ±3cr rejection criterion has been applied to re- 
move high background intervals. After filtering 
for the additional GTIs, the level 2 event files be- 
longing to the same target were co-added. To this 
end, we first improved the absolute astrometry 
(the overall 90% uncertainty circle of Chandra X- 
ray absolute position has a radius of 0.6 arcsec) 
applying corrections for translation, scale and ro- 
tation to the world coordinate system (WCS) of 



Table 1: Chandra observations log for the three 
pointings of the GWS field. Col. (1): observa- 
tion number; Col. (2,3): RA and DEC (J2000); 
Col. (4): Standard Data Processing (SDP) ver- 
sion; Col. (5): Observation format in the Time 
exposure mode; Col. (6): Net exposure time (Sum 
of Good Time Intervals, GTI) (ksec) 



Obs 


RA 


DEC 


SDP 


Mode 


GTI 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


5853 


14:15:22.5 


+52:08:26.4 


7.6.7.1 


VFAINT 


42.4 


5854 






7.6.7.1 


VFAINT 


50.2 


6222 






7.6.7.1 


VFAINT 


35.2 


6223 






7.6.7.1 


VFAINT 


49 


6366 






7.6.7.1 


VFAINT 


15.96 


7187 






7.6.7.1 


VFAINT 


8 


5851 


14:16:24.5 


+52:20:02.59 


7.6.7.1 


VFAINT 


36 


5852 






7.6.4 


VFAINT 


11.95 


6220 






7.6.7.1 


VFAINT 


36 


6221 






7.6.7.1 


VFAINT 


3.8 


6391 






7.6.7.1 


VFAINT 


9.4 


7169 






7.6.4.1 


VFAINT 


18.2 


7181 






7.6.7.1 


VFAINT 


8 


7188 






7.6.4.1 


VFAINT 


4.2 


7236 






7.6.4 


VFAINT 


18.8 


7237 






7.6.4.1 


VFAINT 


17.2 


7238 






7.6.4 


VFAINT 


9.6 


7239 






7.6.4.1 


VFAINT 


16.2 


3305 


14:17:43.0 


+52:28:25.2 


6.8.0 


FAINT 


27.6 


4357 






6.9.0 


FAINT 


79.8 


4365 






6.9.0 


FAINT 


58.2 



each data file, by comparing two sets of source 
lists from the same sky region. Average posi- 
tional accuracy improvement is 6.3% with respect 
to the original astrometry provided by the space- 
craft attitude files. Once we have improved the 
astrometry, the event files have been merged. 

The output merged files have been filtered 
to create event files for several energy bands: 
full (0.5-7keV), soft (0.5-2keV), hard (2-7keV) 
hard2 (2-4.5 keV) and vhard (4-7 keV). Unbinncd 
effective exposure maps at a single energy, rep- 
resentative of each band, were created at 2.5 keV 
(full), IkeV (soft), 4keV (hard), 3keV (hard2) 
and 5.5 keV (vhard). Resulting maps have units of 
time x effective area, i.e. s cm 2 counts photon -1 . 
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2.1.2. Source detection 

We applied the CIAO wavdetect Mexican-Hat 
wavelet source detection 1 program to all bands 
(full, soft, hard, hard2 and vhard) for the three 
fields. Several wavelet scales have been applied: 
1, V2, 2, 2V2, 4, 4V2, 8, 8V2 and 16 pixels. 
Small scales are well suited to the detection of 
small sources while larger scales are appropriate 
for more extended sources. Due to the "blank 
field" nature of the GWS, we do not expect to 
detect sources at a scale larger than 16 pixels (<~ 
8 arcsec). For each applied scale, wtransf orm pro- 
duces a complete set of detections. We set a signif- 
icance threshold of 2 x 10~ 7 . Since this wavdetect 
input parameter measures the number of spurious 
events per pixel, we expect 0.2 fake detections 
per detector of 1024 x 1024 pixels i.e. about 12 
spurious sources in the entire catalog (4 detectors 
x 5 bands x 3 fields). The output wavdetect 
detection files are fits tables containing the posi- 
tion of sources (RA & DEC), count rates in pho- 
tons s _1 cm -2 and ancillary information. The final 
X-ray source catalog has been obtained by cross- 
matching detections in all five bands in order to 
increase the reliability of detections. Best match is 
given by a minimum separation of sources around 
a great circle. We set a maximum search distance 
of 2 arcsec (see section 2.4). We computed the 
hardness ratios as follows: 



HR(A 1 E,A 2 E) 



CR(A!E)-CR (A 2 E) 
CR(A!E) + CR (A 2 E) 



(1) 



where AiE y A 2 E are different energy bands and 
CR(A n E) is the count rate in a given energy 
band. Four hardness ratios have been set in this 
way: HRi = HR(hard, soft), HR^ = HR(hard2, 
soft), HR 2 = HR(vcry hard, hard) and HR' 2 = 
HR(very hard, hard2). 

We also computed the energetic fluxes from count 
rates. To this end, mean photon energies in 
each band have been calculated assuming a power 
law spectrum with T = 1.5 (a reasonable assump- 
tion after inspecting the hardness ratios diagram 
and diagnosis grid depicted in Figure 5), redshift 
z = 0.5, galactic absorption n# = 1.3 x 10 20 cm~ 2 



lr The fact that the PSF of X-ray detectors often has 
Gaussian-like shape motivates our use of the Marr wavelet, 
or "Mexican Hat" function. 



and no intrinsic absorption. We applied the ab- 
sorption coefficients provided by the TRANNM func- 
tion within the package PIMMS v3.9b (Mukai, 
1993). The output catalog has 639 unique X-ray 
emitters (after removing 20 common sources in the 
field overlapping area). From them, 76 objects are 
located outside the field covered by optical data. 
From the total sources, 46% (297) have been de- 
tected with a significance > 4cr, 21% (132) above 
the 3er significance level and 13% (81) down to 
2.5cr level. 

The number of sources detected in each band 
are 490 (full), 465 (soft), 284 (hard), 266 (hard2) 
and 121 (vhard), with median fluxes of 1.81 x 10~ 15 , 
5.57 x 10~ 16 , 2.34 xl0~ 15 , 1.58 x 10" 15 and 
3.75 x 10 -15 erg cm" 2 s _1 , and with limiting fluxes 
at 3cr level of 4.8 x 10~ 16 , 1.1 x 10~ 16 , 2.8 x 10~ 16 , 
1.3xl0~ 16 and 7.3 x 10~ 16 crgcuT 2 s" 1 , respec- 
tively. Median errors in fluxes (for objects above 
the 3a significance level) are 8% in full band, 10% 
in soft band, 12% in hard band, 14% in hard2 band 
and 22% in vhard band. There are 16 sources de- 
tected only in hard band, 18 in hard2 band and 
3 sources are only observed in vhard band. We 
find an overabundance of soft X-ray sources in 
our final catalog. There are 103 sources detected 
only in this band. This effect can be at least par- 
tially explained considering the decrease of tele- 
scope + ACIS effective area from some 600 cm 2 in 
the 1-2 keV range (including the 0.5-2 keV band 
defined as soft), until approximately 200cm 2 at 
~ 6 keV. On the other hand, the total background 
between 0.5 and 7keV is around 4 times higher 
than in the 0.5-2 keV range. The combination of 
these two factors justifies the high Chandra/ ACIS 
source detection efficiency in the soft X-ray band 
and the mentioned overabundance. The first ver- 
sion of this catalog was presented in Sanchez- 
Portal ct al. (2007). 

2.2. Comparison with AEGIS data 

The All-wavelength Extended Groth strip In- 
ternational Survey (AEGIS) collaboration has re- 
cently made public their data sets, including a 
catalog of X-ray sources that comprises the ob- 
servations used within this paper. The complete 
AEGIS X-ray catalog (Laird et al. 2008) consists 
of 1325 sources; from these, 471 are located within 
the area covered by our data (fields EGS-6, 7 and 
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8 according to their nomenclature) . From our 429 
sources detected above 3ct significance level, 416 
coincide with AEGIS sources using a search ra- 
dius of 2 arcsec (388 sources if only unique matches 
are considered). Figure 1 shows a good agreement 
between the source counts used in this work and 
those provided by the AEGIS team: median rela- 
tive deviations of our data with respect to AEGIS 
are 2% (full band), 0.7% (soft), 5.5% (hard) and 
7.8% (vhard). 




200 400 



600 1000 1200 

counts fu n (this work) 




(counts^nfthis workVcountSfonfAEGIS)) x lOO/countSfunfthis work) 



Fig. 1. — Upper panncl: full band net counts from 
this work vs AEGIS net counts. The solid line rep- 
resents the best linear fit, with a slope of 0.99 and 
zero point 5.23 counts. Lower pannel: histogram 
of relative deviations of our full band count data 
with respect to AEGIS. Median deviation is lower 
than 2%. 



2.3. Optical data 

We observed three pointings in the direction of 
the GWS, each of them using the B, V, R, and I 
broadband filters. The total area covered is 0.18 
square degrees. The observations have been car- 



ried out using the Prime Focus Imaging Platform 
(PFIP) at the 4.2m William-Herschel Telescope 
(WHT) in the Roque de los Muchachos Observa- 
tory (La Palma, Canary Islands). This camera has 
two CCD detectors (2k x 4k), with a total field of 
view of 16' x 16' and a pixel size of 0.237 arcsec. 
Several exposures have been performed at each 
pointing, of 600, 800, 900, or 1000 sec, depending 
on the filter, with a dithering of 15 arcsec between 
consecutive exposures in order to allow eliminat- 
ing cosmic rays and to fill the gap between de- 
tectors. Standard reduction procedures have been 
applied. Resulting limiting magnitudes are 25, 25, 
24.5 and 23.5 in B, V, R and I band, respectively. 
A detailed description of data reduction and op- 
tical source detection can be found in Cepa et al. 
(2008). The optical catalog used in this work con- 
tains data for ~ 44000 objects. 

2.4. The catalog of optical counterparts 

In order to build a sample of X-ray emitters 
with optical counterparts, we performed cross- 
matching of the X-ray and optical catalogs using 
the TOPCAT tool (Taylor 2005), applying a crite- 
rion of minimum sky distance along a great circle. 
X-ray sources located inside the field covered by 
optical data and their optical counterparts are rep- 
resented in Figure 2. 

We chose a maximum search radius of 2 arcsec af- 
ter carrying out several tests considering a maxi- 
mum searching distance of 1 to 5 arcsec with a bin 
of 0.5 arcsec (Figure 3). As expected, the num- 
ber of counterparts increases with the maximum 
search radius, but also the number of multiple 
matches within the error circle. Setting a search 
radius of 1 arcsec we minimise the number of mul- 
tiple cross-matches, but at the cost of loosing an 
outstanding number of potential counterparts. On 
the other hand, choosing a search radius greater 
than 2 arcsec produces a large number of multiple 
matches, increasing the uncertainty in counterpart 
identification. Therefore, we set a 2arsec search 
radius as an adequate compromise. This practi- 
cal criterion is consistent with the more rigorous 
analysis performed below. 

We detected 340 optical counterparts of the X- 
ray sources in at least one optical band. Of them, 
322 are detected in the B band, 332 in V, 334 in 
R and 317 in the I band. 

We carried out a quantitative analysis of 
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Fig. 2. — X-ray sources (diamonds) and their opti- 
cal counterparts (filled diamonds) from the cross- 
match of X-ray and optical catalogs 
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Fig. 3. — Number of optical counterparts as a 
function of the maximum search radius used for 
cross-matching. The number of closest sources 
is represented with stars. Crosses represent the 
total number of sources found within the search 
radius. The dashed line marks the adopted sepa- 
ration. The same value was obtained applying the 
statistical methodology of de Ruiter et al. (1977) 

the completeness 1 and reliability 2 of our cross- 
matching procedure by applying the methodology 
of de Ruiter et al. (1977). The separation between 
sources is defined by the dimensionless variable r 
as: 



1 Fraction of all real associations between the X— ray and op- 
tical sources that are indeed classified as identifications. 

2 Fraction of claimed identifications expected to be true 
counterparts. 



(2) 



where Aa and A5 are the position differences be- 
tween the X-ray source and its optical counter- 
part in RA & DEC, and a 2 a and aj are the over- 
all astrometric errors computed as a = [a\ + 
Uopt) 1 ^ 2 - In this work, we assume ax =0.7arcsec 
and a pt =0.3 arcsec. ax was calculated by square 
adding a systematic uncertainity of 0.5 and a cen- 
troiding error of 0.5 arcsec. 
The likchood ratio is defined as: 



LR{r) 



dp(r\id) 1 
dp(r\c) ~ 2A 



e 2 



(2A-1) 



(3) 



where: 



• dp(r|id) is the a priori probability that an 
X-ray source and its optical counterpart 
have a separation between r, and r + dr due 
to astrometrical errors 

• dp(r|c) is the probability that a confusing 
background optical object is found in the 
range r, and r + dr from the X-ray source 
position. 

• A measures the number of confusing ob- 
jects within the error circle and depends on 
the optical objects density (p op t) as: A = 
■na a asp opt . 

An optical source can be considered as a true 
counterpart of an X-ray source if its LR is higher 
than a certain threshold value L. L is obtained 
maximising the sum of completeness and reliabil- 
ity that are, in turn, mathematically defined in 
terms of the probabilities defined above. In our 
case, L = 0.56 and its corresponding variable is 
r = 2.6 which gives us a separation distance of 2.0 
arcsec, a value consistent with the adopted one. 
The completeness of our detection is 98.3% (i.e. 
some 1 object lost due to the use of a particular 
cutoff value L), while reliability is 87.3% (~ 39 
potentially false matchings). 



7 



3. Analysis 



3.1. X ray properties 

We have computed the cumulative number 
count distributions (logN-logS) per deg 2 in the soft 
band and compared it with other surveys in order 
check the reliability of our source detections. We 
have chosen the soft band to be able to compare 
our distribution with that computed by Nandra 
et al. (2005) and those provided by other deep X- 
ray surveys, Chandra Deep Field South (CDFS, 
Giacconi et al. 2002), Chandra Deep Field North 
(CDFN, Alexander et al. 2003), ELAIS-N1 and 
ELAIS-N2 (Manners et al. 2003). All distributions 
are represented in Figure 4. For our sample two 
curves have been computed, one for the complete 
sample of X-ray detections and the other only for 
the X-ray detections that have optical counter- 
parts. The median flux levels are 1.81 x 10~ 15 and 
1.97 Xl0~ ergcm~ 2 s~ 1 , respectively. For the 
objects brighter than ~ 10~ 14 ergs -1 cm~ 2 both 
distributions coincide, while for objects fainter 
than ~ 10~ 14 ergs -1 cm -2 the density of X-ray 
sources with optical counterparts starts to de- 
crease, since the optical sample is incomplete at 
the depth of our X-ray data. In general, at our 
flux limit the number counts of our complete X— 
ray sample are in good agreement with Nandra 
et al. (2005) and other surveys. We have de- 
tected a higher density of sources with flux below 
~ 5 x 10~ 15 erg cm -2 s _1 , most likely due to the 
use of different detection algorithms and more 
sensitive detection threshold. As we already men- 
tioned, our threshold is 2 x 10 -7 , two times higher 
than the one chosen by Nandra et al. (2005) (see 
2.1.2 for more information). 

Three independent energy bands, soft, hard2 
and vhard (see the Section 2.1.1 for ranges), have 
been used to obtain the X-ray color-color diagram 
shown in Figure 5. Hardness ratios have been 
calculated as described in Section 2.1.2. X-ray 
sources are superimposed on a model grid com- 
puted using power law Spectral Energy Distribu- 
tions (SEDs) with different photon indexes, T, af- 
fected by a galactic absorption column density, 
Nff, of 1.3 x 10 20 cm~ 2 and different values of in- 
trinsic absorption; the Chandra ACIS-I response 
matrices have been also included in the computa- 
tion. To obtain the fluxes we have used XSPEC 
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Fig. 4. — Cumulative logN-logS functions for the 
Groth field in the soft band, for all detected X-ray 
sources (thick solid line) and for X-ray sources 
with optical counterpart (thin solid line). Five 
distributions, with different exposure times and 
effective areas, have also been represented as com- 
parison: GWS by Nandra et al. (2005) (dotted 
line), CDFN (thick dash-dot-dash line), CDFS 
(thin dash-dot-dash line), ELAIS-N1 (dash line) 
and ELAIS-N2 (dash-dot-dot-dash line) 

(Arnaud 1996), varying the photon index between 
and 3 with a step size of 0.5, and the intrinsic 
absorption column density between 20 and 24 (in 
logarithmic scale) in steps of 0.5. The majority 
of our sources are located within the model grid 
region, suggesting that their X-ray SED is com- 
patible with the assumed model. A specially high 
concentration of objects can be seen in the region 
of Nh between 10 20 and 5 x 10 21 cm -2 and with 
r~1.8. A group of sources are located outside 
the model grid, showing a soft excess, suggesting 
that they can not be properly described by a single 
power law. For the fitting of these objects proba- 
bly more components should be added. 

3.2. Optical structural parameters 

We obtained a set of structural parameters of 
the optical counterparts of our X-ray sample, in 
order to describe quantitatively their structure 
and morphology. In this section we describe the 
procedure used to derive the parameters and to 
produce a morphological classification. 

We used SExtractor v2.2.2 (Bertin & Arnouts 
1996) to detect and extract the sources. Details 
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Fig. 5. — X-ray color-color diagram based on 
2-4.5/0.5-2keV (HRi) vs. 4-7/2-4.5keV (HR' 2 ) 
hardness ratios. X-ray sources (crosses), from the 
complete X-ray catalog, are represented and su- 
perposed on an absorbed power law model grid. 
Photon index is varying between and 3 in steps of 
0.5 (dashed lines from up to down) and a log(N#) 
between 20 and 24 in steps of 0.5 (solid lines from 
left to right). Median error bars are represented 
in the upper left corner of the diagram. 

of the source extraction process are provided in 
Cepa et al. (2008). The version used (v2.2.2) was 
modified by B. Holwerda to measure the Abra- 
ham concentration index (Abraham et al. 1994). 
This is defined as the ratio between the integrated 
flux within certain radius defined by the normal- 
ized radius a and the total flux. We used the 
values obtained at a = 0.3 in our analysis in or- 
der to compare our results with those of Abraham 
et al. (1996). Another relevant output parame- 
ters for our study are the CLASS-STAR, which 
permits to roughly classify objects as compact or 
extended (for more information see Section 3.2.1), 
and the background that will be used later. 

In addittion, we used the Galaxy IMage 2D 
software package (GIM2D, Simard 1998; Simard 
et al. 2002) versions 2.2.1 and 3.1, to perform the 
galaxy structural parameter decomposition (see 
Figure 6). This software considers two compo- 
nents: an spheroid, represented by a Sersic pro- 
file 2 , and a disk represented by an exponential 
profile. Using the standard steps in DAOPHOT 



2 While any rigourous analysis of the structural components 
of an AGN host should include an additional component, 
namely a pointlikc nucleus, the small size of most of our 



(Stetson 1987) we obtained the PSF images re- 
quired for GIM2D execution. The initial pa- 
rameters (including background estimation) and 
setup for GIM2D were determined using the im- 
ages and catalogs produced by SExtractor (Bertin 
& Arnouts 1996). We obtained from the output 
provided by GIM2D the following parameters 3 : 
bulgc-to-total ratio (B /T) , residual parameter (R) 
(Schade et al. 1995, 1996), asymmetry (A) in- 
dex (Abraham et al. 1994, 1996), the Sersic index 
(n) and the total luminosity {L tot ). The resid- 
ual parameter and the asymmetry index quantify 
the galaxy irregularity, by measuring the deviation 
from the simple, symmetric model (e.g. spheroid 
+ disk), using as input the residual and origi- 
nal image, respectively. GIM2D also calculates 
for certain output parameters their lower and up- 
per limits at 99% confidence level. In order to 
estimate the accuracy of the best fit parameters 
we have looked at the deviation, computed as the 
maximum between the best value and the upper 
or lower limits (see Figure 7). In this respect, the 
B/T parameter and Sersic index have revealed in- 
accurate when applied to faint and small objects 
(< 300 pixel, while the seeing disk (FWHM) area 
is about 30 pixels) as shown in Figure 7. Therefore 
we have decided to do not use these parameters for 
the morphological classification. 

In practice, we observed that the code either 
fails or produces unreliable output parameters in 
a number of cases: (1) when the source R magni- 
tude is less than 18 or larger than 24, (2) when the 
isophotal area is very small (usually less then 90 
pixels), (3) when the central object is surrounded 
by close neighbouring sources or nearby bright 
companions, and (4) when the source is close to 
the frame boundaries. 



objects compared with the seeing disk prevents from imple- 
menting such a detailed approach. Rather, we have simply 
assumed that the sum of the nuclear and bulge components 
can be roughly represented by a single Sersic profile. 
3 According to the available documentation (Simard 1998; 
Simard ct al. 2002), GIM2D calculates the Abraham con- 
centration at four normalized radii a = 0.1, 0.2, 0.3 and 
0.4. (Abraham ct al. 1994, 1996). Comparing the re- 
sults with the ones obtained with SExtractor, a systematic 
shift of ~0.2 is found towards larger values in GIM2D fig- 
ures. After performing several tests we concluded that the 
GIM2D code performs its computations using normalized 
area rather than normalized radius against what is stated 
in the documentation. 
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Fig. 6. — Four sample images of the GIM2D 
decomposition procedure. First column shows 
subimages (obtained from the scientific image) of 
objects to be fitted; the second column repre- 
sents the mask subimages, third column shows 
the subimages obtained from the modelling and 
the fourth column shows the residual images, ob- 
tained from the original image after subtracting 
the model 
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Fig. 7. — Evolution of the median error of the 
bulgc-to-total ratio B/T (crosses), Sersic index n 
(stars) and total flux (filled circles) L^ ^ as func- 
tions of the isophotal area. Deviations are evalu- 
ated as explained in Section 3.2. 

3.2.1. Morphological classification 

The selection of compact objects (QSOs, faint 
host AGNs or active stars) has been accomplished 
by using the SExtractor CLASS_STAR output pa- 
rameter. The computed value depends on the 
seeing, the peak intensity of the source and its 
isophotal area. The parameter ranges between 



(extended object like a galaxy) and 1 (pointlike 
object). A source has been deemed as compact if 
its CLASS_STAR parameter is larger than 0.9. 

To determine the morphology of extended ob- 
jects (i.e. not selected as compact according to 
the criterion above), we analyzed four parame- 
ters: Abraham concentration index (C), Abraham 
asymmetry index (A), B/T ratio and residual pa- 
rameter (R). We carried out the analysis in three 
different bands (B, V, and R) in order to determine 
the most reliable parameters, and the best band 
for morphological classification. We excluded I 
band here due to fringing. 

We first classified each object using the set of 
parameters mentioned above, separately for each 
parameter in each band, considering the results 
from previous works as limits between different 
morphological types. We used the results of Im 
et al. (2002) and Simard (1998) to obtain prelimi- 
nary morphology with the B/T parameter, Schade 
et al. (1995) results for the residual parameter and 
Abraham et al. (1996) for the classification with 
the concentration and asymmetry index. In prac- 
tice, we found that, in many cases, different pa- 
rameter sets lead to a different classification of 
a given object even in the same band. In order 
to have an independent evaluation criterion, we 
also carried out a visual classification of the ob- 
jects. Despite its subjectiveness, this procedure 
can yield highest-quality results when applied to 
bright and extended objects, thus being a conve- 
nient method to test the reliability of the struc- 
tural parameters that will be in turn applied to 
the analysis of fainter objects. The visual classifi- 
cation has been accomplished by inspection of the 
radial profile, 2D surface plot and isophotal coun- 
tours diagram of each source (see Figure 8). For 
the dimmest and smallest objects of our sample 
(galaxies with the magnitude R > 24 and isopho- 
tal area < 100 pixels) the profiles and isophotal 
contours are too noisy to get an acceptable result 
by means of a simple visual inspection. 

From this comparison we conclude that for our 
sample, characterised by faint objects (typically R 
> 22) with small isophotal area (< 300 pixel, while 
the seeing FWHM disk area is about 30 pixels) the 
most reliable parameter is the concentration index 
C, combined either with the assymmetry index or 
the residual parameter. We therefore based our 
morphological classification in the concentration 
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index combined with both asymmetry index and 
residual parameter. The R band has been used, 
since, on the one hand, the largest part of optical 
counterparts have been found in this band, and on 
the other, the achieved S/N is generally larger. B, 
V, and I bands have been used to perform the mor- 
phological classification of objects not detected in 
the R band, or detected in this band but with very 
low S/N ratio. 




Fig. 8. — Sample of different Hubble type galax- 
ies and their contour diagrams; Up: Left: Group 
I (E/SO), Center: group II (Early type spiral), 
Right: group III (Spiral); Down: Left, group III 
(Spiral), Center: group IV (Irr) , Right: compact 
object (C>0.7 or CLASS.STAR > 0.9). All ob- 
jects are located in the center of the images. 

Based on the classification using visual inspec- 
tion of the brightest sample objects, and upon 
comparison with previous works (e.g. Abraham et 
al. 1996, see references above) in the present pa- 
per we created five categories or groups associating 
ranges of values of C with Hubble types: 

0: C > 0.7 or CLASS_STAR > 0.9 -» Compact obj. 
I: 0.45<C<0.7^E, E/SO and SO 
II: 0.3 < C < 0.45 -> SO/SOa-Sa 
III: 0.15<C<0.3^Sab-Scd 
IV: C<0.15^Sdm-Irr 

The first group is populated by elliptical and 
lenticular galaxies, whose profiles do not show any 
trace of disk or spiral arms (see Figure 8), the sec- 
ond group includes very early spiral type galax- 
ies which show signs of spiral arms. The third 



group comprises spiral galaxies while the fourth 
includes late type spirals and irregular galaxies, 
that show almost no regular form in their profiles. 
All galaxies, even those in the first group (ellip- 
tical and lenticular) have a concentration index 
C < 0.7, while all objects with C > 0.7 have been 
classified by SExtractor as compact ones, with the 
CLASS.STAR parameter > 0.9. 

In Figure 9 we depict the Abraham asymmetry 
vs. concentration index diagram segregated ac- 
cording to the morphology determined using the 
visual classification for the brightest objects, and 
the structural parameters (in most of the cases C 
combined with A) for the fainter ones. Compar- 
ing our Figure 9 with that shown in Abraham et 
al. (1996), it can be seen that in our case there 
is a small shift of the concentration index towards 
smaller values. This can be explained consider- 
ing the influence of the atmospheric seeing con- 
ditions. SExtractor does not apply any correction 
on the computed parameters for this effect, and as 
the seeing induced PSF increases, both concentra- 
tion and asymmetry indexes tend to decrease, as 
confirmed by our simulations carried out by con- 
volving simple galaxy models with seeing PSFs of 
different widths. Since the Abraham et al. (1996) 
diagram has been obtained from HST data, these 
differences are expected. 

The morphological analysis described above has 
been applied to the complete sample of 340 X-ray 
objects with optical counterpart (see Figure 10). 
Of those, 333 were classified using the R band, 4 
sources using the V band and finally 3 objects by 
means of B band data. The largest group (143 
sources) is composed by compact objects (QSOs, 
faint host AGNs, stars) followed by group III (73 
galaxies). 57 objects are found either in group II 
or III, having outstanding bulge component a signs 
of a formed disk and spiral arms. Only 33 galax- 
ies are found to belong to the IV group. Finally, 
we detected 7 closed pairs, whose possible phys- 
ical relation will be confirmed when photometric 
redshifts are computed (paper III, in preparation). 
Due to the high noise level, or to the presence of 
other objects close to the observed ones, the clas- 
sification uncertainty for some of the objects is 
quite high, and it was not possible to determine 
the morphological group: 3 galaxies could be ei- 
ther in group II or III and other 3 either in group 
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Fig. 9. — Morphological classification diagram 
combining the Abraham asymmetry index (A) 
vs. Abraham concentration index (C). Diamonds 
present E and SO galaxies, triangles SO/SOa-Sa 
galaxies, squares young spiral galaxies (Sab-Sbcd) 
and crosses late type galaxies (Sdm-Irr). Com- 
pact objects are represented by star symbols. Dot- 
ted lines separate the locii of elliptical, spiral and 
peculiar galaxies (from right to left) obtained by 
Abraham et al. (1996). There is a small shift of 
our concentration and asymmetry index towards 
smaller values due to the influence of the atmo- 
spheric seeing conditions. 

Ill or IV. For 21 sources it was not possible to 
determine the morphology since they are at the 
detection limit. 

4. Relationship between X-ray and opti- 
cal structural properties 

4.1. Nuclear classification 

The simplest way to classify the nuclear type is 
based on the X-ray-to-optical ratio (X/O). The 
typical value of X/O flux ratios for X-ray se- 
lected AGN (both type 1 or BLAGN and type 
2 or NLAGN) is in the range between 0.1 and 
10 (e.g. Fiore et al. 2003). For optically selected 
type 1 AGNs the typical value is X/O ~ 1 (Alexan- 
der et al. 2001; Fiore et al. 2003). At high X/O 
flux ratios (well above 10) we can find BLAGN 
and NLAGN as well as high-z high-luminosity ob- 
scured AGN (type 2 QSOs), high-z clusters of 
galaxies and extreme BL Lac objects. Finally, the 
X/O < 0.1 region is typically populated by coronal 
emitting stars, normal galaxies (both early type 



Fig. 10. — The histogram represents the final mor- 
phological classification of the GWS X-ray emit- 
ters with optical counterparts. Different values of 
the Abraham concentration index correspond to 
different groups I-IV (see text). A group of possi- 
ble interacting objects has been detected; a group 
of sources (marked with '?') remained unidenti- 
fied due to their very low S/N ratio. 

and star-forming) and nearby heavily absorbed 
(Compton thick) AGN. 

We have applied a simple criterion for perform- 
ing a coarse nuclear-type classification of our sam- 
ple objects, based on diagnostic diagrams relating 
X-ray-to-optical ratios (X/O) to hardness ratios 
(HR) . This criterion is based on the results of Delia 
Ceca et al. (2004) on the XMM-Newton Bright 
Serendipitous Survey; they found that, when plot- 
ting the HR(2 - 45fceVy0.5 - 2keV) (equivalent 
to our HR'j) against X/O, computed as the ra- 
tio of the observed X-ray flux in the 0.5-4.5keV 
(soft+hard2) energy range and the optical R-band 
flux, most (85%) of the BLAGN identified by 
means of optical spectroscopy are tightly packed 
in a small rectangular region of the diagram, while 
NLAGN (type 2) tend to populate a wide area of 
the diagram (towards harder values of the hard- 
ness ratio). We applied the same approach to our 
sample, combining the HR^ hardness ratio with 
the X/O ratio, where the optical flux Fr has been 
derived from the Petrosian magnitudes in the R 
band. The resulting diagram, that also includes 
host morphology information, is shown in Figure 
11, where the dashed-line box indicates the locus 
of the BLAGN region. We set the right edge of 
this box, HR(2 - 4.5fceVy0.5 - 2fceK)«-0.35 as 
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a coarse boundary separating 'unobscured' AGN 
(BLAGN) and 'obscured' AGN (NLAGN). On the 
other hand, we consider the small box depicted in 
the diagram as the "highest probability" region 
for finding a BLAGN (Delia Ceca et al. 2004). 
We find that a large fraction of our objects, 63% 
4 , fall inside the region of BLAGNs, while 51% 
of the sample is placed within the small BLAGN 
box defined by Delia Ceca et al. (2004). Regard- 
ing morphology, we do not find evidence of any 
relationship between nuclear and morphological 
types, since we did not obtain any clear separation 
of nuclear species based on morphological type. 
Nevertheless, we can see that the majority of ob- 
jects identified as compact (65%) are placed in the 
BLAGN region. It is also interesting to highlight 
the relationship between the X/O ratio and the 
host galaxy morphology, with early-type galaxies 
having generally lower values for X/O ratio and 
an opposite behavior of late-type objects. Finally, 
we find that only about 7% of our sample objects 
are placed in the lower region of the diagram (i.e. 
that typical of stars/normal galaxies/Compton- 
thick AGN). 

We also represented the X/O ratio with respect 
to HRi, yielding a basically identical distribution 
as that described above, and with respect to HR 2 
and HR' 2 , observing that in these two cases is 
harder to get a clear separation between differ- 
ent nuclear/morphological types; this behavior is 
expected since hardness ratios computed from the 
hard and vhard bands are less sensitive to absorp- 
tion than those involving the soft band, which is 
one of the main criteria for the separation between 
broad and narrow line AGN. 

4.2. X ray and optical properties 

We combined the X/O ratio with optical col- 
ors (B-R and B-I) for the different morphologi- 
cal groups and nuclear types, as shown in Fig- 
ure 12. We observe that compact and late-type 
objects tend to show bluer colors, while early- 
type galaxies tend to present redder colors and 
also lower X/O values, as expected. On the other 
hand, there is a mild tendency for BLAGN to show 



4 A11 fractions given here are relative to the total number 
of X— ray objects having optical counterpart and detected 
in both soft and hard2 energy bands such as HR(2 — 
4.5fceV/0.5 — 2keV) can be computed 



bluer colors than those found in NLAGN, as al- 
ready observed in local universe samples (e.g. Yee 
1983; Sanchez-Portal et al. 2004). The B-R av- 
erage value for BLAGNs is 1.0 and for NLAGNs 
1.4, with a — 0.5 and 0.6 respectively. Using the 
Kolmogorov-Smirnov test, the two distributions 
are statistically different with a significance level 
of 99%. 

When representing the X/O ratio vs. the Abra- 
ham concentration index C we find a strong anti- 
correlation between these two parameters (Figure 
13), obtaining a correlation coefficient of ~-0.5. 
Using Spearman and Kendall (S-K) statistics, we 
obtained a correlation significance level >99%. 
The first order polynomial function has been fit- 
ted with a slope of -2.03, logX/O-intercept of -0.94 
and a standard deviation of <~0.5. This result is 
consistent with that referred above and depicted 
in Figure 11, that is the tendency of early- type 
galaxies to have lower than average X/O values 
and the opposite behavior in late-type ones. As 
seen in Figure 13, all morphological groups, in- 
cluding compact objects, tend to follow the same 
relation. For compact objects we obtained simi- 
lar parameters using S-K statistics: a correlation 
coefficient of ^0.5 and a correlation significance 
level >99%. A least squares fit gives to the anti- 
correlation a slope of -2.4 and a logX/O-intercept 
of -0.95. Comparing different nuclear types, we 
found that both BLAGN and NLAGN follow the 
same distribution (see Figure 13). On the other 
hand, we did not find any clear correlation be- 
tween HR(2 - 4.5/ceU/0.5 - 2keV) and the con- 
centration index (Figure 14), suggesting that the 
observed anticorrelation is not related to obscura- 
tion. 

The physical origin of the observed anticorrela- 
tion is not well established. First of all, it is not 
easy to disentangle possible bias effects that super- 
imposed to variable seeing conditions. It has been 
pointed out the existence of redshift-dependent bi- 
ases in both A and C, due to the reduction of ap- 
parent size and surface brightness with respect to 
the sky background with increasing z, along with 
poorer sampling and lower S/N (Brinchmann et 
al. 1998; Conselice et al. 2003); on the other hand, 
there is a "bandpass shift" with respect to the rest 
wavelength (Bershady et al. 2000). Even if these 
effects combine to produce some uncertainty in the 
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Fig. 11. — Relationship between X/O ratio and HR(2 — 4.5keV/0.5 — 2keV), for the different morphological 
types. Solid lines separate BLAGN and NLAGN regions and the area with X/O < 0.1, where stars, normal 
galaxies and Compton thick AGNs can be found (Fiore et al. 2003). Dashed line box presents the limits 
obtained by Delia Ceca et al. (2004) where ~ 85% of their spectroscopically identified BLAGN have been 
found. In our case, inside of this box can be found ~ 51% of all objects from the diagram. 
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Fig. 12. — Log(X/0) flux ratio vs. B-R color for different morphology (left) and nuclear types (right). 



morphological classification of galaxies, generally 
in the sense of shifting objects to apparently later 
Hubble types (i.e. with lower C values), A and 
C seems to be reproducible out to z = 3, within a 
reasonable scatter (Conselice et al. 2003). More- 
over, there can be a truly evolutionary effect in 
the AGN host galaxies, indicating that galaxies at 
higher redshift arc intrinsically less concentrated 
(Grogin et al. 2005). 

At this stage we cannot apply any correction 
which depends on the redshift, although this is 
one of the objectives for the future work. How- 
ever, we have performed a preliminary evaluation 



by representing the concentration index C versus 
the spectroscopic redshifts z gathered from the 
DEEP2 catalog for 84 objects from our sample. 
While for non compact objects we obtained an an- 
ticorrelation between these two parameters, C and 
z are found uncorrelated for compact objects, in 
clear contrast with the observed behavior in Fig- 
ure 13. 

So, although the biases described above might 
contribute partially to it, the observed anticorre- 
lation seems to also be a result of a true physical 
effect. In order to attempt to explain it, we have 
first to understand the physical connections of the 
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Fig. 13. — Relationship between X/O and the Abraham concentration index C for different morphological 
(left) and nuclear types (right). A clear anticorrelation is observed. 



1.0 



-0.5 



-1.0 



- * Compact objects 




- o E/S0 




- a Early type spirals 


\ o 


Spiral galaxies 


- + Sdm/lrr 








+ 


□ 






— □ 


_ * a 

* * 




X 




4* * * ° 


+ 


□ H„ ^ * „ ^ 


+ 


* O 

x m. „* * * * * x 








if S ** 




* * A^K» * 
m O * 




A * 



-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 -0.0 

logC 



Fig. 14. — Relationship between the hardness ratio HR(2 — 4.5keV/0.5 — 2keV), and the concentration index 
C for the different morphological types. 



two observables C and X/O. The concentration 
index (C) is tightly linearly correlated to the cen- 
tral velocity dispersion (Graham et al. 2001a,b), 
which is in turn related to the nuclear black hole 
(BH) mass (e.g. Magorrian et al. 1998). There- 
fore, the concentration and the nuclear BH mass 
(Mbh) are correlated. The physical connections 
of the X-ray to optical flux ratio are not so ev- 
ident. The X/O ratio measures the X-ray flux 
(in the 0.5-4.5keV band), normalized with the R- 
band optical flux of the whole galaxy (nucleus, 
bulge and disk) . It should be noted that no band- 
pass, or K-correction is performed; therefore if 
the sample spans a large redshift range the in- 
terpretation of the quantity becomes more prob- 



lematic. If the nuclear luminosity is large com- 
pared with that of the host galaxy, X/O ratio can 
be thought as a measure of the X-ray to optical 
spectral index, aox 5 - This parameter has been 
found to be strongly anticorrelated with the UV 
luminosity ^ 2500 ^i without a significant correla- 
tion with redshift (Steffen et al. 2006). On the 
other hand Bian (2005) found a strong correla- 
tion between the hard/soft X-ray spectral index, 
ax and the Eddington ratio in a sample of 41 
BLAGN and narrow-line Seyfert 1 (NLS1) galax- 
ies observed with ASCA. Other works have con- 



5 This relation is usually defined as aox 
0.3838 l° S (h keV /f,) 
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firmed this correlation (Grupe 2004; Shemmer et 
al. 2006). The relationship between ax and aox 
has not been thoroughly studied in large sam- 
ples of AGN, but a linear correlation has been 
found between both spectral indexes in a sample 
comprising 22 out of 23 quasars in the complete 
the Palomar-Green X-ray sample with z < 0.4 and 
Mb < -23 (Shang et al. 2007). If this last correla- 
tion holds for our sample, X/O could be tracing 
the Eddington ratio in the large nuclear luminosity 
limit. In addition, if the host galaxy luminosity is 
large when compared with the nuclear luminosity 
(Lbuige+disk > L nuc i eus ), the X/O ratio can be 
thought lower bound of the X-ray-to-bulge 
luminosity ratio, that is in turn a (weak) measure 
of the AGN Eddington ratio L/L Edd cx L/M B h, 
assuming that the X-ray luminosity represents 
the nuclear luminosity and the bulge luminosity 
is proportional to the bulge mass (and therefore 
to M M ). 

Hence, we can guess a (loose) correlation be- 
tween X/O and the energy production efficiency 
of the AGN. Under these assumptions, the C 
vs. X/O relation traces the correlation between 
the Eddington ratio and the nuclear BH mass. 
The obtained result could therefore suggest that 
early-type galaxies, having poor matter supply to 
feed the activity, have lower Eddington rates than 
those of late-type, with larger reserves of the gas 
for AGN feeding. 

This suggested approach is consistent with 
the results obtained by Ballo et al. (2006, 2007) 
in a sample of X-ray selected AGNs at z<l. 
They observed that a large fraction of low- 
luminosity AGNs are fed by black holes with 
a mass M B h > 3 x 10 6 MQ, and Eddington rate 
< 1, while a small ( < 10%) fraction of AGNs 
host small SMBH (M B h < 10 6 M©) and Edding- 
ton rate w 1, which corresponds to a more effi- 
cient accretion rate. They stated that their results 
strongly suggest that most of the low-luminosity 
X-ray selected AGNs (at least up to z < 0.8) are 
powered by rather massive black holes, experienc- 
ing a low-level accretion in a gas-poor environ- 
ment. 

Our approach is also consistent with the re- 
sults of Wu & Liu (2004), who studied the black 
hole masses and Eddington rates of a sample of 
135 double-peaked broad line AGNs, observed in 
two surveys: SDSS and a survey of BLAGN ra- 



dio emitters. They obtained that if the separa- 
tion between the line peaks (which is correlated 
to the SMBH mass) decrease, the Eddington rate 
increases. 

Kawakatu et al. (2007) found an anticorrelation 
between the mass of a SMBH and the luminosity 
ratio of infrared to active galactic nuclei Edding- 
ton luminosity, L/^/L^d, for a sample of ultralu- 
minous infrared galaxies with type 1 Scyfcrt nuclei 
(type 1 ULIRGs) and nearby QSOs, which also 
could be consistent with our approach. The an- 
ticorrelation is interpreted as a link between the 
mass of a SMBH and the rate of mass accretion 
onto a SMBH, normalized by the AGN Eddington 
rate, which indicates that the growth of the black 
hole is determinated by the external mass sup- 
ply process, and not the AGN Eddington-limited 
mechanism, changing its mass accretion rate from 
super-Eddington to sub-Eddington. 

5. The clustering of X-Ray emitters and 
extended optical counterparts 

The study of large-scale environmental condi- 
tions of AGN host galaxies are useful not only 
for understanding the differential properties of the 
AGN phenomena, but to trace the formation and 
evolution of such galaxies and the structures where 
they reside. Despite the relevance of X-ray se- 
lected AGNs, their clustering properties still re- 
main quite unknown (Basilakos et al. 2005, and 
references therein). In this sense, and taking into 
account the characteristics of the X-ray source 
catalog presented in this work, we have calculated 
the 2-point angular correlation function (2p-ACF) 
for significant subsamples. 

The 2p-ACF, u)(9), gives the excess of probabil- 
ity, with respect to a random homogeneous distri- 
bution, of finding two sources in the solid angles 
dfii and cif^ separated by an angle 9, and it is 
defined as 



dP=M 2 dfl 1 dQ 2 



(4) 



where N is the mean number density of cataloged 
sources. To measure the 2p-ACF we have follow 
the same procedure as done in Cepa et al. (2008), 
adopting the estimator of Landy & Szalay (1993), 
which can be written in the form 



DD - 2DR + RR 
RR 



(5) 
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where DD — N ss (9i) is the fraction of possible 
source-source pairs counted in i-bins over the an- 



gular range studied and DR = 



(N g -1) 



2N r 



N gr (9 t 



RR 



N,(N,-1) 
N r {N r -l) 



N rr (6i) 



are the normalized 
counts of source-random and random-random 
pairs in these bins, respectively. To make the ran- 
dom catalogs we placed (200 x N source ) random 
points by following the X-ray source angular dis- 
tribution. Howewer, the source density is affected 
by several instrumental biases such as detector 
gaps and sensivity variations due to survey re- 
gions with different exposure times. To account 
for these effects, we made normalized sensitivity 
maps from the effective exposure time and source 
distributions for each band used in this section 
(soft and hard2). The raw maps were modelled 
using a mimimum-curvature algorithm. Next, the 
uniform density random catalogs were convolved 
with the corresponding sensitivity map to obtain 
differential number counts in the surveyed area. 
In the case of the HR(2 - A.5keV/0.5 - 2keV) 
catalog, we adopt the sensitivity map used for the 
hard2 source one. The smoothed footprint of each 
sensitivity map can be seen in Figure 15. 

We measure u>(9i) in i— logarithmic bins of 
width A(log 9) = 0.2 in the angular scales 0.24 < 
6 < 24.0 arcmin. The lower limit in angular sepa- 
ration is according to the minimum detection scale 
defined in the Section 2.1.2. The formal error as- 
sociated with the 2p-ACF measured is the Poisso- 



nian estimate of, = 



i+h)(gi) 



DD . The 2p-ACF binned 
measures were fitted using the A LU (9 1 ~' < — C) power 
law. We adopted a fixed canonical slope 7=1.8 
(the so called "comoving clustering model") and 
C is defined as 



C = 



E Nrr(0) 9~ { 

E N rr (9) 



(6) 



corresponding to the numerical estimation of the 
"integral constraint" (Peebles 1980). As in Miyaji 
et al. (2007), we use the normalization A as the fit- 
ting parameter rather than the correlation length 
9 C = A 1 ^ 7-1 ), since the former gives better con- 
vergence of the fit. 

From the whole catalog we have selected the 
soft (0.5-2 keV) and hard2 (2-4.5 kcV) sources 
(given their similar depths in the survey), the X— 
ray objects with HR(2 - 4.5keV/0.5 - 2keV)^0 
and the fraction of X-ray sources with optical 
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Fig. 15. — Smoothed sensitivity maps in the soft 
(0.5-2 keV) and hard2 (2-4.5 keV) bands, used 
to convolve the uniform density random catalogs. 
The sensitivity dimension is normalized to the 
maximum exposure time in each band. Contours 
are linearly binned and lighter regions correspond 
to higher sensitivities. 



counterparts. Table 2 summarizes the basic results 
of the correlation analysis for these subsamples. 
In all cases, we have obtained significant positive 
clustering signal. Regarding the soft band sam- 
ple, the correlation length 9 C derived from A w is 
consistent with the obtained by Basilakos et al. 
(2005) (XMM- Newton LSS) and it is - 2 o larger 
than the values obtained by Gandhi et al. (2006), 
Carrera et al. (2007) and Ueda et al. (2008) . In the 
specific case of the hard2 band, we only found in 
the literature the 2p-ACF estimations of Miyaji et 
al. (2007), valid for X-ray sources detected in the 
XMM-Newton observations of the COSMOS field 
and corresponding to their MED band: our fitted 
correlation length is ~4.3 times larger, although 
our sample is a factor 2 deeper. 

In Figure 16 are represented our 2p-ACF mea- 
surements with the best-fit A^ 9~ - 8 power-law 
in soft and hard2 bands (panels a and b, respec- 
tively). It is noticeable that the clustering ampli- 
tude obtained in both bands looks similar. Thus, 
in the case of the Groth field, we would not expect 
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Table 2: Clustering properties of the X-ray selected sources. N source is the number of data points used in 
the 2p-ACF estimation and A u is the best- fit value in the separation range given by # m ; n — 6* max , assuming 
a power-law u>(0) = A u #( 1_ t) =_0 - 8 . The goodness-of-fit statistics (x 2 ) and the fitted 2p-ACF at 1 arcmin 
[w(l')] are given in the final columns. The errors represent la uncertainties. 



Subsample 


N 

iv source 


^min ^max 


A u (1-7 = -0.8) 


x 2 


u>{\>) 






arcmin 


xl0~ 3 degrees 






soft (0.5-2 keV) 


465 


0.24-24.0 


6.99 ±0.68 


0.030 


0.15 ±0.02 


hard2 (2-4.5 keV) 


266 


0.24-24.0 


9.19 ± 1.28 


0.107 


0.20 ±0.03 


HR(2-4.5fceV/0.5-2fceV) 


190 


0.24-24.0 


9.68 ± 1.30 


0.109 


0.21 ±0.03 


Optical counterparts (I<23) 


104 


0.60-24.0 


19.99 ±3.32 


0.096 


0.42 ±0.07 



redshift or spatial clustering differential effects in 
the soft and hard2 band samples. Interestingly, 
an absence of correlation in differently defined 
hard bands has been reported by Gandhi et al. 

(2006) and Puccetti et al. (2006); nevertheless, 
Yang et al. (2003) and Basilakos et al. (2004) not 
only yield a positive correlation amplitude, but 
the hard band signal is several times that of the 
soft band. Other results with similar correlation 
lengths in both bands are given by Carrera et al. 

(2007) , Miyaji et al. (2007) and Yang et al. (2006). 

In the panel c of Figure 16 is represented the 2p- 
ACF estimations fo a HR(2 - 4.5keV/0.5 - 2keV) 
selected sample. The clustering amplitude is 
slightly larger than the soft or hard2 band ones. 
Moreover, the whole HR(2 - 4.5fceV/0.5 - 2keV) 
subsample was decomposed using the bound- 
ary referred in Sect. 4.1 to separate type 1 (or 
BLAGN) from type 2 (or NLAGN) populations 
and a positive, very similar clustering signal from 
both subsamples emerged from the 2p-ACF es- 
timations, despite their goodness-of-fit statistics 
are worse than the tabulated values in a factor 
~ 1 because of the small size of the subsamples. 

Finally, we measured the 2p-ACF of the op- 
tical counterparts of X-ray sources in order to 
compare the clustering behavior with those that 
come from the analysis of galaxy populations in 
the same field. For a reliable comparison is nec- 
essary to diskard the compact objects (i.e. with 
CLASS_STAR>0.9) from the optical counterparts 
catalog, as well as all objects fainter than a lim- 
iting magnitude. We have selected the peak of 
the apparent brightness distribution in the I-band 
as such limiting magnitude. Thus, X-ray coun- 
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Fig. 16. — 2p-ACF estimates of the X-ray sources 
detected in the soft band (panel a), hard2 band 
(panel b) and HR as defined in cq. 1 (panel c). 
The continuous line represents the A u #( 1_ t)=~ - 8 
power law fit to the data. 

terparts with I>23 were ruled out, remaining 
104 extended objects, which can be regarded as 
AGN host candidates of the corresponding X-ray 
sources. 

In spite of this skinny sample, it is possible to 
obtain some interesting lines of thought about the 
X-ray selected AGN clustering. It is suggested 
that the typical environments of AGNs are no dif- 
ferent to those of inactive galaxies in general, at 
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redshift ranges above z^0.4 (e.g. Grogin et al. 
2003; Waskett et al. 2005). Nevertheless, Coil 
ct al. (2004) and Gilli et al. (2005) have argued 
that the AGN are preferentially hosted by early- 
type galaxies at a typical redshift of <~ 1 . In this 
way, Basilakos et al. (2004) found that clustering 
length values of hard X-ray objects are compara- 
ble with those of Extremely Red Objects (EROs) 
and luminous radio sources. More recently, Geor- 
gakakis et al. (2007) demonstrated that X-ray- 
selected AGNs at z<~ 1, in the AEGIS field, avoids 
undcrdense regions at 99.89% of significance. 

In Figure 17 we have depicted the 2p-ACF 
of the optical counterparts selected as described 
above. The 2p-ACFs of the Groth field galaxy 
sample with I<23 only and with V-I>3 (I<24) 
from Cepa et al. (2008) are overplotted. At least 
in the range 0.01 < 9 <0.1 degrees, the cluster- 
ing signal produced by the optical counterparts of 
our X-ray sources is l-2cr stronger than the corre- 
sponding to the I< 23 galaxy population. In terms 
of fitted correlation length, the 9 C (X-ray) (1-7=- 
0.8) is ~ 30 times larger than that of the Groth 
field galaxies, but only ~ 2.3 times larger than the 
maximum correlation length of the V-I>3 galax- 
ies in the same field. In the other hand, using 
the optical data from the counterparts catalog, we 
found that the median V-I observed color is 1.86 
and most of 65% of the optical counterparts have 
an observed color V-I>2. This means that a sig- 
nificant population of our optical counterparts is a 
composition of red and very red galaxies, whereas 
its clustering signal tends to reproduce that asso- 
ciated with such galaxies, which drive the ^10x 
clustering excess found by Cepa et al. (2008). At 
once, more than a 55% of the optical counterparts 
belong to the morphological categories I and II, as 
defined in the Section 3.2.1. 

6. Summary and Conclusions 

We processed and analyzed public, 200 ksec 
Chandra / ACIS observations of three fields com- 
prising the original GWS gathered from the Chan- 
dra Data Archive, combined with our broadband 
survey carried out with the 4.2m WHT at La 
Palma. Our X-ray catalog contains 639 unique 
X-ray emitters, with a limit flux at 3cr level of 
4.8 x 10 -16 erg cm" 2 s _1 and a median flux of 1.81 
x 10~ 15 ergcm _2 s _1 in the full (0.5-7 keV) band. 
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Fig. 17. — 2p-ACF estimation for the opti- 
cal extended (CLASS.STAR<0.9) counterparts 
with I magnitude brighter than 23 (solid cir- 
cles). The continuous thick line represents the 
A w 0(!-7)=-o.8 p 0W e r l aw fit to the data. The 
dashed line represents the 2p-ACF estimation for 
galaxies with I<23 in the Groth field, whereas the 
continous thin line it does for the V-I>3 color se- 
lected galaxies with I < 24, both extracted from 
Cepa et al. (2008). The shaded regions represent 
fitting errors of lcr. 

Cross-matching the X-ray catalog with our opti- 
cal broadband catalog, we have found 340 X-ray 
emitters with optical counterparts. The complete- 
ness and the reliability of cross-matching proce- 
dure are 97% and 88.9%, respectively. At our 
flux limit, the number counts of our complete X- 
ray sample are in good agreement with Nandra 
et al. (2005) and other comparison surveys men- 
tioned in Section 3.1. X-ray SEDs of the mayor- 
ity of our sources are compatible with the power 
law model, affected with the galactic absorption 
column density and the intrinsic absorption. We 
have performed a morphological classification ob- 
taining and analyzing different structural param- 
eters, including the bulge- to- total ratio, residual 
parameter, Sersic index, asymmetry and concen- 
tration indexes. We also performed the classifica- 
tion based on visual inspection. To obtain the nu- 
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clear classification, we applied a criterium based 
on the correlation between X/O flux ratio and 
HR1' hardness ratio. 63% of our sources, with 
HRl'^0, have been classified as BLAGNs. The 
electronic table summarizes the X/O flux ratios, 
HR1' values and morphological classification of all 
340 X-ray emitters with optical counterparts. We 
have investigated the correlations between X-ray 
and broadband optical structural parameters that 
provide with information about the host galaxy 
morphology and populations, finding that: 

1. For a sample where most of the objects are 
faint (R > 22) and with small isophotal area (< 300 
pixels), bulge-to-total flux ratio has revealed inac- 
curate to perform the morphological classification. 
Concentration index combined with the asymme- 
try index have been found to provide the best 
separation criterion for such a sample, different 
morphological types have been observed to dis- 
tribute according to the following ranges of con- 
centration index: group (compact objects) — > 
C > 0.7 (or CLASS_STAR > 0.9), group I (E, E/S0 
and SO) -» 0.45<C<0.7, group II (SO/SOa-Sa) 
-» 0.3<C<0.45 (0.47), group III (Sab-Scd) 
0.15 < C < 0.3 and group IV (Sdm-Irr) -» C < 0.15 

2. No clear separation has been found between 
the morphology and the nuclear type, but most of 
our objects classified as compact are placed in the 
region of BLAGNs. 

3. We can confirm the tendency of compact and 
late type objects to show bluer colors, while early 
type galaxies tend to present redder colors. On the 
other hand, we find a mild tendency of BLAGNs to 
show bluer colors when compared with NLAGNs. 

4. We obtained an anticorrelation between the 
X-to-optical flux ratio and the concentration in- 
dex (which parametrizes the morphology), with 
significance level higher than 99%, showing that 
the early-type galaxies tend to have lower and the 
late- type galaxies higher than average X/ O values. 
Objects classified as compact show this anticor- 
relation as well, with the same significance level. 
The anticorrelation might suggest that early-type 
galaxies, having poor matter supply to feed the ac- 
tivity, have lower accretion rates than those of late 
type, with larger reserves of gas for AGN feeding. 

Finally, we presented the first results on the 
two-point angular correlation function of the X- 
ray selected sources in the Groth field from Chan- 



dra /ACIS observations, segregated in four sub- 
samples: soft, hard2, HR(2_4.5fc e v/o.5— 2fceV) an d 
optical extended counterparts. We obtain the fol- 
lowing results and conclusions: 

1. All subsamples reveal a significant positive 
clustering signal in the ^0.5 to 6 arcmin separa- 
tion regime. Canonical power-law (7=1.8) fits to 
the angular correlation functions have been made, 
obtaining correlation lengths between ~ 7 and 27 
arcsec, depending on the subsample. 

2. We have neither found substantial differ- 
ences between the clustering of soft and hard2 
selected X-ray sources, nor when comparing the 
populations of type 1 and type 2 AGN (i.e. 'unob- 
scured' and 'obscured' AGN, respectively), sepa- 
rated using the HR(2-4.5£:eVyO.5-2fceV0 w - 0.35 
boundary. 

3. The clustering analysis of the optical coun- 
terparts of our X-ray emitters provides a correla- 
tion length significantly larger than the one found 
for the whole galaxy population with the same lim- 
iting magnitude, and same field, but similar to 
that of strongly clustered populations like red and 
very red (bulge-dominated) galaxies. A clustering 
signal > 1(7 was observed even at 6 ~ 30 arc- 
sec, which corresponds to 160 h^ 1 kpc, assuming 
a mean redshift of 0.85 for Chandra extragalactic 
sources (Hasan 2007). This comoving distance is, 
for instance, comparable to the separation of two 
mergers in the t«2 Gyr stage of their large-scale 
evolution, according to the model of Mayer et al. 
(2007) on rapid formation of SMBH. Hence, our 
results suggest that the environment plays an im- 
portant (but not unique) role as possible triggering 
mechanism of AGN phenomena. This fact can be 
perfectly compatible with the interpretation given 
to the anticorrelation between X/O ratio and the 
concentration index. The former may be associ- 
ated to the genesis of the AGNs, while the anti- 
correlation could be that to their maintenance and 
evolution. 
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